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Abstract 



Consider transportation of one distribution of mass onto another, chosen to optimize 
the total expected cost, where cost per unit mass transported from x to y is given by 
C^ . a smooth function c{x,y). If the source density f^{x) is bounded away from zero and 

infinity in an open region U' C R", and the target density f~{y) is bounded away 
from zero and infinity on its support V C R", which is strongly c-convex with respect 
to U', and the transportation cost c satisfies the (A3)w condition of Trudinger and 
^ . Wang |51| , we deduce local Holder continuity and injectivity of the optimal map inside 

■"ijij- ' U' (so that the associated potential u belongs to C'jq" ([/')). Here the exponent a > 

depends only on the dimension and the bounds on the densities, but not on c. Our result 
provides a crucial step in the low/interior regularity setting: in a sequel [T7|, we use it 
to establish regularity of optimal maps with respect to the Riemannian distance squared 

C^ ' on arbitrary products of spheres. Three key tools are introduced in the present paper. 

Namely, we first find a transformation that under (A3)w makes c-convex functions level- 
set convex (as was also obtained independently from us by Liu |40|1. We then derive new 
Alexandrov type estimates for the level-set convex c-convex functions, and a topological 
lemma showing optimal maps do not mix interior with boundary. This topological lemma, 
which does not require (AS)^, is needed by Figalli and Loeper [20] to conclude continuity 
^ ' of optimal maps in two dimensions. In higher dimensions, if the densities /^ are Holder 

continuous, our result permits continuous differentiability of the map inside U' (in fact, 
^loc regularity of the associated potential) to be deduced from the work of Liu, Trudinger 
and Wang pT]. 
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1 Introduction 

Given probability densities < /^ G L^(R") with respect to Lebesgue measure ^" on R"-, 
and a cost function c : R" x R" i — > [0, +oo], Monge's transportation problem is to find a 
map G : R" i — > R" pushing d//"*" = f^dJ^"^ forward to d/i" = f~d^''^ which minimizes the 
expected transportation cost |46| 



inf I c{x,G{x))dfi-^{x), (1.1) 



G 

where G^fi^ = /.f~ means ^~[y] = fi'^[G~^{Y)] for each Borel Y C R". 

In this context it is interesting to know when a map attaining this infimum exists; sufficient 
conditions for this were found by Gangbo [25] and by Levin [39], extending work of a number 
of authors described in |26| |54| . One may also ask when G will be smooth, in which case it 
must satisfy the prescribed Jacobian equation | deti^G(a;)| = f^{x)/f~{G{x)), which turns 
out to reduce to a degenerate elliptic partial differential equation of Monge- Ampere type for 
a scalar potential u satisfying Du{x) = —Dxc{x,G{x)). Sufficient conditions for this were 
discovered by Ma, Trudinger and Wang |45| and Trudinger and Wang |51| |52| . after results 
for the special case c{x,y) = \x — yp/2 had been worked out by Brenier [1], Delanoe |12| . 
Caffarelli [6] [5] [7] [8] [9], and Urbas |53| . and for the cost c(x, y) = — log |x — y| and measures 
supported on the unit sphere by Wang |56| . 



If the ratio f^{x)/f~{y) — although bounded away from zero and infinity — is not 
continuous, the map G will not generally be differentiable, though one may still hope for it 
to be continuous. This question is not merely of technical interest, since discontinuities in / 
arise unavoidably in applications such as partial transport problems |10| [3] |15| |16| . Such 
results were established for the classical cost c{x,y) = \x — yp/2 by Caffarelli [5] [7] [8], for 
its restriction to the product of the boundaries of two strongly convex sets by Gangbo and 
McCann |27| . and for more general costs satisfying the strong regularity hypothesis (A3) 
of Ma, Trudinger and Wang |45) — which excludes the cost c{x,y) = \x — y|^/2 — by 
Loeper |42) ; see also |34| |40| |52) . Under the weaker hypothesis (AS)^ of Trudinger and 
Wang |51j . which includes the cost c{x,y) = \x — yp/2 (and whose necessity for regularity 
was shown by Loeper |42| . see also |23|). such a result remained absent from the literature; 
the aim of this paper is to fill this gap, see Theorem 12.11 below. 

A number of interesting cost functions do satisfy hypothesis (A3)w, and have applica- 
tions in economics |18| and statistics |48) . Examples include the Euclidean distance squared 
between two convex graphs over two sufficiently convex sets in R" |45| . the simple harmonic 
oscillator action |38| . and the Riemannian distance squared on the following spaces: the round 
sphere |43| and perturbations thereof |13| |21| |22| , multiple products of round spheres (and 
their Riemannian submersion quotients, including products of complex projective spaces) [36], 
and products of perturbed 2-dimensional spheres [ll]. As remarked in |35j , for graphs which 
fail to be strongly convex and all Riemannian product geometries, the stronger condition 
(A3) necessarily fails. In a sequel, we apply the techniques developed here to deduce regu- 
larity of optimal maps for the multiple products of round spheres [17]. Moreover, Theorem 12.1 
allows one to apply the higher interior regularity results established by Liu, Trudinger and 
Wang |41[ . ensuring in particular that the transport map is C°°-smooth if f~^ and /~ are. 



Most of the regularity results quoted above derive from one of two approaches. The 
continuity method, used by Delanoe, Urbas, Ma, Trudinger and Wang, is a time-honored 
technique for solving nonlinear equations. Here one perturbs a manifestly soluble problem 
(such as I detZ?Go(a;)| = f~^{x)/fo{Go{x)) with /o = f~^, so that Go{x) = x) to the problem 
of interest (| det-DGi(x)| = f~^{x)/ fi{Gi{x)), /i = /~) along a family {ft}t designed to 
ensure the set of i € [0, 1] for which it is soluble is both open and closed. Openness follows 
from linearization and non-degenerate ellipticity using an implicit function theorem. For the 
non-degenerate ellipticity and closedness, it is required to establish estimates on the size of 
derivatives of the solutions (assuming such solutions exist) which depend only on information 
known a priori about the data (c, /i). In this way one obtains smoothness of the solution 
y = Gi{x) from the same argument which shows Gi to exist. 

The alternative approach relies on first knowing existence and uniqueness of a Borel map 
which solves the problem in great generality, and then deducing continuity or smoothness 
by close examination of this map after imposing additional conditions on the data (c, / ). 
Although precursors can be traced back to Alexandrov [2] , in the present context this method 
was largely developed and refined by Caffarelli [5] [7] [H], who used convexity of u crucially 
to localize the map G{x) = Du{x) and renormalize its behaviour near a point [x,G{x)) of 
interest in the borderline case c{x,y) = —{x,y). For non-borderline (A3) costs, simpler 



estimates suffice to deduce continuity of G, as in |27| |11) |42) |52) : in this case Loeper was 
actually able to deduce an explicit bound a = (An — 1)~^ on the Holder exponent of G when 
n > 1. This bound was recently improved to its sharp value a = (2n— 1)~^ by Liu |40| . using 
a key observation discovered independently from us (see Section U] and Theorem 14. 3p ; both 
Loeper and Liu also obtained explicit exponents a = a{n,p) for /"*" G L^ with p > n |32] or 
p > (n + l)/2 |40| and l//~ G -L°°. For the classical case c{x,y) = —{x,y), explicit bounds 
were found by Forzani and Maldonado |24) . depending on log f-(l G L'^ [H] |55| . 

Our proof introduces at least three significant new tools. Its starting point is that condi- 
tion (A3)w allows one to add a null Lagrangian term to the cost function and exploit diffeo- 
morphism (i.e. gauge) invariance to choose coordinates, which depend on the point of interest, 
to transform the c-convex functions to level-set convex functions (see Theorem I4.3p ; this ob- 
servation was made also by Liu [40J, independently from us. Next we establish Alexandrov 
type estimates (Theorems 16.21 and 16. lip for c-convex functions whose level sets are convex, 
extending classical estimates for convex functions. These rely on quantitative new aspects 
of the geometry of convex sets that we derive elsewhere, and which may have independent 
interest |19j . The resulting Alexandrov type estimates enable us to exploit Caffarelli's ap- 
proach [5] [7] [8] more systematically than Liu |40j was able to do, to prove continuity and 
injectivity (see Theorems 18. II and 18.2^ . Once such results are established, the same estimates 
permit us to exploit Forzani and Maldonado's |23] approach to extend the engulfing prop- 
erty of Gutierrez and Huang |29) to (A3)w c-convex functions (see Theorem 19. 3p . improving 
mere continuity of optimal maps to Holder continuity (Theorem 19.51 and its Corollary 19. 6p . 
Along the way, we also have to overcome another serious difficulty, namely the fact that 
the domain of the cost function (where it is smooth and satisfies appropriate cross-curvature 
conditions) may not be the whole of R". (This situation arises, for example, when optimal 
transportation occurs between domains in Riemannian manifolds for the distance squared 
cost or similar type.) This is handled by using Theorem 15. H where it is first established that 
optimal transport does not send interior points to boundary points, and vice versa, under 
the strong c-convexity hypothesis (B2)s described in the next section. Theorem 15.11 does 
not require (A3)w, however. Let us point out that, in two dimensions, there is an alter- 
nate approach to establishing continuity of optimal maps; it was carried out by Figalli and 
Loeper |20| following Alexandrov's strategy [2], but their result relies on our Theorem 15.11 

2 Main result 

Let us begin by formulating the relevant hypothesis on the cost function c{x, y) in a slightly 
different format than Ma, Trudinger and Wang |45| |51| . We denote their condition (A3)w 
as (B3) below, and their stronger condition (A3) as (B3)s. Throughout the paper, Dy will 
denote the derivative with respect to the variable y, and iterated subscripts as in D"^ denote 
iterated derivatives. For each {x,y) ^U xV we define the following conditions: 

(BO) [/ C R" and F C R" are open and bounded and c G C^(I7 x F); 



(Bl) (bi-twist) — n^ f~\ ( ^^^ diffeomorphisms onto their ranges; 

y (z V I >■ UxC\x,y) J 

Uy := -Dyc{U,y) \ _^ ___ ,^^^^ „„„. 



(B2) (bi-convex) ^}^ ' t^ )~ \r\ f are convex subsets of R' 
Vx:= -Dxc{x,V) J 

(B3) (-(A3)w) for every curve t G [-1,1] ^ {Dyc[x{t),ym,D.,c{x{Q),y{t))) G R^" 

which is an afiinely parameterized hne segment, 



9^ 

cross(3.(o),y(o))[x'(0),y'(0)] := 



provided 



c(x(s),y(i))>0 (2.1) 

(s,t)={0,0) 



d' 



dsdt 



c{x{s),y{t)) = 0. (2.2) 

(s,t)={0,0) 



From time to time we may strengthen these hypotheses by writing either: 

(B2)s if the convex domains Uy and Vx in (B2) are strongly convex; 

(B3)s (=(A3)) if, in the condition (B3), the inequahty (|2.ip is strict; 

(B4) if, in the condition (B3), (I2.ip holds even in the absence of the extra assumption (12.2 



Here a convex set Q C R" is said to be strongly convex if there exists a radius R < +oo 
(depending only on Q), such that each boundary point x £ dQ can be touched from outside 
by a sphere of radius R enclosing Q; i.e. Q C Br{x — Rhg^x)) where nQ{x) is an outer 
unit normal to a hyperplane supporting Q at x. When Q is smooth, this means all principal 
curvatures of its boundary are bounded below by 1/R. 

Hereafter U denotes the closure of U , int U denotes its interior, diam U its diameter, and 
for any measure //+ > on [/, we use the term support and the notation spt ^~^ C U to refer 
to the smallest closed set carrying the full mass of fi~^. 

Condition (A3)w (=(B3)) was used by Trudinger and Wang to show smoothness of 
optimal maps in the Monge transportation problem (jl.ip when the densities are smooth. 
Necessity of Trudinger and Wang's condition for continuity was shown by Loeper |42| . who 
noted its covariance (as did |35) |49| ) and some relations to curvature. Their condition 
relaxes the hypothesis (A3) (=(B3)s) proposed earlier with Ma |45| . In |35| . Kim and 
McCann showed that the expressions (12. ip and (12. 2p correspond to pseudo-Riemannian sec- 
tional curvature conditions induced by the cost conU xV, highlighting their invariance under 
reparametrization of either [/ or F by diffeomorphism; see |35| Lemma 4.5]; see also [37] as 
well as [30] . for further investigation of the pseudo-Riemannian aspects of optimal maps. 
The convexity of Uy required in (B2) is called c-convexity of U with respect to y by Ma, 
Trudinger and Wang (or strong c-convexity if (B2)s holds); they call curves x{s) G U, for 
which s G [0,1] I — 7- —Dyc{x{s),y) is a line segment, c-segments with respect to y. Similarly, 
V is said to be strongly c*-convex with respect to x — or with respect to U when it holds for 
all X £ U — and the curve y(t) from (|2.ip is said to be a c*-segment with respect to x. Such 



curves correspond to geodesies (x{t),y) and {x,y{t)) in the geometry of Kim and McCann. 
Here and throughout, line segments are always presumed to be affinely parameterized. 
We are now in a position to summarize our main result: 

Theorem 2.1 (Interior Holder continuity and injectivity of optimal maps). Let c £ C'^{U x 
V) satisfy (B0)-(B3) and (B2)s. Fix probability densities /+ E L^{U) and /" e i^(^) 
with {f^ / f~) GL°°([/xy) and set d^ji^ := f^d^^ . (Note that spt fi'^ may not be c-convex.) 
If the ratio {f~ /f~^) G L°°{U' x V) for some open set U' C U (U' is not necessarily c-convex), 
then the minimum U.l\} is attained by a map G : U i — > V whose restriction to U' is locally 
Holder continuous and one-to-one. Moreover, the Holder exponent depends only on n and 
l|log(/+//^)||Loo((//xV)- 

Proof. As recalled below in Section[3](or see e.g. |54| ) it is well-known by Kantorovich duality 
that the optimal joint measure 7 G r(^''~, /i~) from p.ip vanishes outside the c-subdifferential 
(j3.3p of a potential u = u^ ^ satisfying the c-convexity hypothesis (j3.2p . and that the map 
G : U I — > V which we seek is uniquely recovered from this potential using the diffeomorphism 
(Bl) to solve (jS.Sp . Thus the Holder continuity claimed in Theorem 12.11 is equivalent to 

^ e ci^:{u'). 

Since /U^ do not charge the boundaries of U (or of V), Lemma lS-lT e) shows the c-Monge- 
Ampere measure defined in (j3.6p has density satisfying l^'^ul < \\f~^ /f~\\L°°(UxV) o^^ ^ ^^^ 

ll/"//+IIZk^'xV) ^ l^'^l ^ Wf^/f'h^iU'xV) on U'. Thus u G GJ^^^iU') according to 
Theorem 19.51 Injectivity of G follows from Theorem 18. H and the fact that the graph of G is 
contained in the set d'^u C U x V of (j3.3p . The dependency of the Holder exponent a only 
on n and || log(/^//~)||/,oo(f//x\/) follows by Corollary 19.61 D 

Note that in case /+ G Gc{U) is continuous and compactly supported, choosing U' = 
U^ = {f~^ > e} for all e > 0, yields local Holder continuity and injectivity of the optimal 
map y = G{x) throughout Uq. 

Theorem 12.11 allows to extend the higher interior regularity results established by Liu, 
Trudinger and Wang in [H], originally given for (A3) costs, to the weaker and degenerate 
case (A3)w, see |411 Remark 4.1]. Note that these interior regularity results can be applied to 
manifolds, after getting suitable stay-away-from-the-cut-locus results: this is accomplished 
for multiple products of round spheres in [T7], to yield the first regularity result that we 
know for optimal maps on Riemannian manifolds which are not flat, yet have some vanishing 
sectional curvatures. 

Let us also point out that different strengthenings of the (A3)w condition have been 
considered in |45| |35) |36| |44| |21) |22| . In particular, one stronger condition is the so-called 
non-negative cross-curvature, which here is denoted here by (B4). Although not strictly 
needed for this paper, under the (B4) condition the cost exponential coordinates introduced 
in Section [3] allow to deduce stronger conclusions with almost no extra effort, and these 
results play a crucial role in the proof of the regularity of optimal maps on multiple products 
of spheres |17] . For this reason, we prefer to include here some of the conclusions that one 
can deduce when (A3)w is replaced by (B4). 



3 Background, notation, and preliminaries 

Kantorovich discerned |32| |33| that Monge's problem (jl.ip could be attacked by studying 
the linear programming problem 



mill / c{x,y)dj{x,y). (3.1) 

Here T{fi^,fi~) consists of the joint probability measures on U x V C R" x R" having // 
for marginals. According to the duality theorem from linear programming, the optimizing 
measures 7 vanish outside the zero set of u{x) + v{y) + c(x, y) > for some pair of functions 
(u, v) = {v'^, u^ ) satisfying 

v'^ix) := sup -c(x, y) - v{y), u"* (y) := sup -c(x, y) - u{x); (3.2) 

these arise as optimizers of the dual program. This zero set is called the c-subdifferential of 
u, and denoted by 



5^= 



{ix,y)GUxV\ u{x) + n^* (y) + c{x, 2/) = o} ; (3.3) 



we also write d^u{x) := {y \ {x,y) G d'^u}, and d'^ u'^ {y) := {x \ {x,y) £ d'^u}, and 
d^u{X) := L)x£xd^u{x) for X C R". Formula (13. 2p defines a generalized Legendre-Fenchel 
transform called the c-transform; any function satisfying u = u^ ^ := {u^ Y is said to be 
c-convex, which reduces to ordinary convexity in the case of the cost c{x,y) = —{x,y). In 
that case d^u reduces to the ordinary subdifferential du of the convex function u, but more 
generally we define 

du := {{x,p) £ U X R" I u{x) > u{x) + {p,x — x) + o{\x — x\) as x — )■ x}, (3.4) 

du(x) := {p I {x,p) G du}, and du{X) := VJx&xdu{x). Assuming c G C^(C/ x F) (which 
is the case if (BO) holds), any c-convex function u = u'^ '^ will be semi-convex, meaning its 
Hessian admits a bound from below D^u > — ||c||(72 in the distributional sense; equivalently, 
u(x) + ||c||(72|xp/2 is convex on each ball in U |26| . In particular, u will be twice-differentiable 
^"-a.e. on U in the sense of Alexandrov. 

As in |25| |39| |45| . hypothesis (Bl) shows the map G : doni Du 1 — > V is uniquely defined 
on the set dom Du C C/ of differentiability for u by 

Dxc{i,G{x)) = -Du{i). (3.5) 

The graph of G, so-defined, lies in d'^u. The task at hand is to show (local) Holder continuity 
and injectivity of G — the former being equivalent to u G Gi^"{U) — by studying the relation 
d^uCU X V. 

To this end, we define a Borel measure \d'^u\ on R" associated to u by 

la^KAT) :=^"(a"n(X)) (3.6) 



for each Borel set X C R"; it will be called the c-Monge- Ampere measure of u. (Similarly, 
we define \du\.) We use the notation \d'^u\ > A on [/' as a shorthand to indicate IS'^tiKX) > 
A^"(X) for each X C C/'; similarly, \d''u\ < A indicates Id'^uKX) < K^''{X). As the next 
lemma shows, uniform bounds above and below on the marginal densities of a probability 
measure 7 vanishing outside d'^u imply similar bounds on \d'^u\. 

Lemma 3.1 (Properties of c-Monge-Ampere measures). Let c satisfy (BO)-(Bl), while u 

and Uj^ denote c-convex functions for each A; € N. Fix x £ X and constants A, A > 0. 

(a) Then d^u{U) C V and \d^u\ is a Borel measure of total mass ^"(y) on U . 

(h) If Uk — ?• Uoo uniformly, then Uoo is c-convex and \d'^Uk\ -^ l^^^^ool weakly-* in the duality 

against continuous functions on U x V . 

(c) Ifuk {x) = for all k, then the functions u^ converge uniformly if and only if the measures 
\d'^Uk\ converge weakly-*. 

(d) If \d''u\ <AonU, then la^'u'^*! > 1/A on V. 

(e) If a probability measure 7 > vanishes outside d^u dU xV , and has marginal densities 
f^, then f~^>XonU'cU and f^<AonV imply |3'^ii| > A/A on U' , whereas f^ < A 
on U' and f" > \ onV imply \d^u\ < A/A on U' . 

Proof, (a) The fact d^u{U) C V is an immediate consequence of definition (|3.3p . Since 
c (z C^iJJ xV), the c-transform v = u'^ : V 1 — t- R defined by ()3.2p can be extended to a 
Lipschitz function on a neighbourhood of V, hence Rademacher's theorem asserts dom Dv is a 
set of full Lebesgue measure in V. Use (Bl) to define the unique solution F : doniDu 1 — > U 
to 

Dyc{F{y),y) = -Dv{y). 

As in |25| |39| . the vanishing of u{x) + v{y) + c{x,y) > implies d'^ v{y) = {F{y)}, at least 
for all points y £ douiDv where V has Lebesgue density greater than one half. For Borel 
X C R", this shows d'^u{X) differs from the Borel set F~^{X) n F by a _Sf" negligible subset 
of V, whence \d'^u\ = F^(^^"'[yj so claim (a) of the lemma is established. 

(b) Let ||ufc — UcxjW 1^00 /jj) — )• 0. It is not hard to deduce c-convexity of iioo, as in e.g. |18| . 
Define v^ = u^* and F^ on domDvk C V as above, so that I^^^Ufcl = Ffc#(^'^[y). Moreover, 
Vk -^ ^00 in L°°{V), where Voo is the c*-dual to «oo- The uniform semiconvexity of Vk (i.e. 
convexity of Vk{y) + ^||c||(72|2/p) ensures pointwise convergence of Dvk — )• Dvoo ^"-a.e. on 
V. From Dyc{Fk{y),y) = —Dvk{y) we deduce Fk — )• -Fqo .i?"-a.e. on V. This is enough 
to conclude jS'^iifel ^ |5'^Mfc|, by testing the convergence against continuous functions and 
applying Lebesgue's dominated convergence theorem. 

(c) To prove the converse, suppose Uk is a sequence of c-convex functions which vanish 
at X and |5'^iifc| ^ fioo weakly-*. Since the Uk have Lipschitz constants dominated by \\c\\(ji 
and U is compact, any subsequence of the Uk admits a convergent further subsequence by the 
Ascoli-Arzela Theorem. A priori, the limit Uoo might depend on the subsequences, but (b) 
guarantees jS'^tiool = l^-oo, after which [421 Proposition 4.1] identifies Uoo uniquely in terms of 
A*^ = Aioo and fi~ = ^""[y, up to an additive constant; this arbitrary additive constant is 
fixed by the condition Uaoix) = 0. Thus the whole sequence Uk converges uniformly. 



(e) Now assume a finite measure 7 > vanishes outside d'^u and has marginal densities 
/^. Then the second marginal dfi~ := f~d^"' of 7 is absolutely continuous with respect 
to Lebesgue and 7 vanishes outside the graph of F : V 1 — > U, whence 7 = (F x id)^iJ,~ 
by e.g. [1] Lemma 2.1]. (Here id denotes the identity map, restricted to the domain domDv 
of definition of F.) Recalling that \d'^u\ = -F#(^"[y) (see the proof of (a) above), for any 
Borel X C U' we have 

X\d'u\{X) = \^^{F-\X)) < [ f-{y)d^^{y) = [ f+{x)d^^{x) < A^^{X) 

Jf-^{x) Jx 

whenever \ < f~ and /+ < A. We can also reverse the last four inequalities and interchange 
A with A to establish claim (e) of the lemma. 

(d) The last point remaining follows from (e) by taking 7 = (F x id)^^'". Indeed an 
upper bound A on \d^u\ = F^^"' throughout U and lower bound 1 on ^" translate into a 
lower bound 1/A on \d^ u^ |, since the reflection 7* defined by 'y*{Y x X) := ^{X x Y) for 
each X xY G U xV vanishes outside d^ u'^ and has second marginal absolutely continuous 
with respect to Lebesgue by the hypothesis \d'^u\ < A. D 

Remark 3.2 (Monge- Ampere type equation). Differentiating (j3.5p formally with respect to 
X and recalling | detDG(x)| = f'^{x)/f~{G{x)) yields the Monge- Ampere type equation 

det[DlMx) + DlxCJx, G{x))] ^ f+{x) 

|detD2^c(x,G(x))| /-(G(x)) ^-^ 

on U, where G{x) is given as a function of x and Du{x) by (13. Sp . Degenerate ellipticity 
follows from the fact that y = G{x) produces equality in u{x) + u'^ (y) + c{x, y) > 0. A 
condition under which c-convex weak-* solutions are known to exist is given by 

[_f+{x)d^^{x)= [_f-{y)d^^{y). 
Ju Jv 

The boundary condition d'^u{U) C V which then guarantees Du to be uniquely determined 
/+-a.e. is built into our definition of c-convexity of u. In fact, [42l Proposition 4.1] shows u 
to be uniquely determined up to additive constant if either /+ > or /^ > .=5f"-a.e. on its 
connected domain, U or V. 

A key result we shall exploit several times is a maximum principle first deduced from 
Trudinger and Wang's work |51| by Loeper; see [42^ Theorem 3.2]. A simple and direct proof, 
and also an extension can be found in |35^ Theorem 4.10], where the principle was also called 
'double- mountain above sliding-mountain' (DASM). Other proofs and extensions appear 
in [52] EOl [Ml EH [21]: 



Theorem 3.3 (Loeper's maximum principle 'DASM'). Assume (B0)-(B3) andfixx,x E 
U. If t £ [0,1] I — > —Dxc{x,y{t)) is a line segment then f{t) := —c{x,y{t)) + c{x,y{t)) < 
max{/(0),/(l)}/ora//tG [0,1]. 



It is through this theorem and the next that hypothesis (B3) and the non-negative cross- 
curvature hypothesis (B4) enter cruciahy. Among the many coroharies Loeper deduced from 
this result, we shall need two. Proved in |42| Theorem 3.1 and Proposition 4.4] (alternately 
|351 Theorem 3.1] and |34| A. 10]), they include the c-convexity of the so-called contact set 
(meaning the c*-subdifferential at a point), and a local to global principle. 

Corollary 3.4. Assume (B0)-(B3) and fix {x,y) £ U xV. If u is c-convex then d'^u{x) 
is c* -convex with respect to x £ U , i.e. —Dxc{x,d^u{x)) forms a convex subset of TiJJ . 
Furthermore, any local minimum of the map x £ U i — > u{x) + c{x, y) is a global minimum. 

As shown in |36| Corollary 2.11], the strengthening (B4) of hypothesis (B3) improves 
the conclusion of Loeper's maximum principle. This improvement asserts that the altitude 
f{t) at each point of the evolving landscape then accelerates as a function oit£ [0, 1]: 

Theorem 3.5 (Time-convex DASM). Assume (B0)-(B4) and fix x,x £U. If t £ 

[0, 1] I — > —Dxc{x,y{t)) is a line segment then the function t £ [0, 1] i — > f{t) := —c{x,y{t))-\- 
c{x,y{t)) is convex. 

Remark 3.6. Since all assumptions (B0)-(B4) on the cost are symmetric in x and y, all 
the results above still hold when the roles of x and y are exchanged. 

c-Monge- Ampere equation 

Fix A, A > and an open domain U^ C U, and let u be a c-convex solution of the c-Monge- 
Ampere equation 

A^" < |9^u| < i^" in U^ C U, , . 

\d'=u\< A^" in U. ^ ^^ 

Note that throughout this paper, we require c-convexity (B2) not of U but only of U. We 
sometimes abbreviate (j3.8p by writing \d^u\ £ [A, 1/A]. In the following sections we will prove 
interior Holder differentiability of u on U , that is u £ Ci^^{U ); see Theorems 18.21 and 19.51 

Convex sets 

We close by recalling two nontrivial results for convex sets. These will be essential in Section [6] 
and later on. The first one is due to Fritz John [Mj: 

Lemma 3.7 (John's lemma). For a compact convex set Q C R" with nonempty interior, 
there exists an affine transformation L : R" — ?• R" such that Bi C L^^{Q) C Bn. 

The above result can be restated by saying that any compact convex set Q with nonempty 
interior contains an ellipsoid E, whose dilation nE by factor n with respect to its center 
contains Q: 

E CQ CnE. (3.9) 

The following is our main result from [19]; it enters crucially in the proof of Theorem 16. Ill 
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Theorem 3.8 (Convex bodies and supporting hyperplanes). Let Q C R" be a well-centered 
convex body, meaning that (j3.9p holds for some ellipsoid E centered at the origin. Fix < 
s < sq < 1. To each y £ (1 — s)dQ corresponds at least one line i through the origin and 
hyperplane 11 supporting Q such that: 11 is orthogonal to i and 

dist(2/, n) < c{n, so)s^/2""' diam(£ n Q). (3.10) 

Here, c{n, sq) is a constant depending only on n and sq, namely c{n, sq) = n^''^{n—-^)( - — fL/2n ' 

4 Choosing coordinates which "level-set convexify" c-convex 
functions 

Recall that c G C^(C/ x y) is a cost function satisfying (B1)-(B3) on a pair of bounded 
domains U and V which are strongly c-convex with respect to each other (B2)s. In the cur- 
rent section, we introduce an important transformation (mixing dependent and independent 
variables) for the cost c{x,y) and potential u{x), which plays a crucial role in the subsequent 
analysis. This change of variables and its most relevant properties are encapsulated in the 
following Definition 14.11 and Theorem 14.31 

Definition 4.1 (Cost-exponential coordinates and apparent properties). Given c G C^{U x 
y) strongly twisted (BO)-(Bl), we refer to the coordinates {q,p) £ Uy x Vx defined by 

q = q{x) = -Dyc{x,y), p = p{y) = -Dxc{x,y), (4.1) 

as the cost exponential coordinates from y £ V and x £ U respectively. We denote the 
inverse diffeomorphisms by x : Uy C T^V i — > U and y : Vx C T~U i — > V ; they satisfy 

q = -Dyc{x{q),y), p = -Dxc{x, y{p)). (4.2) 

The cost c{q,y) = c{x{q),y) — c{x{q),y) is called the modified cost at y. A subset of U 
or function thereon is said to appear from y to have property A, if it has property A when 
expressed in the coordinates q £ Uy. 

Remark 4.2. Identifying the cotangent vector 0©g with the tangent vector q* (BO to U xV 
using the pseudo-metric of Kim and McCann |35) shows x{q) to be the projection to U of the 
pseudo-Riemannian exponential map exp^^ v)il* ® ^)' similarly y(p) is the projection to V of 
expr^y\{0 ©p*). Also, x{q) =: c*-expy q and y{p) =: c-exp^p in the notation of Loeper [32]. 

In the sequel, whenever we use the expression c{q, •) or u{q), we refer to the modified 
cost function defined above and level-set convex potential defined below. Since properties 
(B0)-(B4) (and (B2)s) were shown to be tensorial in nature (i.e. coordinate independent) 
in [35] [42], the modified cost c inherits these properties from the original cost c with one 
exception: (|4.2p defines a C^ diffeomorphism q £ Uy i — > x{q) £ U, so the cost c £ C^{Uy x V) 
may not be C^ smooth. However, its definition reveals that we may still differentiate c four 
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times as long as no more than three of the four derivatives fall on the variable q, and it leads 
to the same geometrical structure (pseudo-Riemannian curvatures, including (j2.ip ) as the 
original cost c since the metric tensor and symplectic form defined in |35) involve only mixed 
derivatives D^ c, and therefore remain C^ functions of the coordinates {q, y) € Uy x V. 
We also use 

(3t = /3t{UxV):= \\iD%c)^'U^^uxV) (4.3) 

^t = 7t{UxV):=\\det{D%c)^'U..^uxV) (4.4) 

to denote the bi-Lipschitz constants /3^ of the coordinate changes (|4.ip and the Jacobian 
bounds 7^ for the same transformation. Notice 7^7(7 > 1 for any cost satisfying (Bl), 
and equality holds whenever the cost function c{x, y) is quadratic. So the parameter 7^7^ 
crudely quantifies the departure from the quadratic case. The inequality /3^(3^ > 1 is much 
more rigid, equality implying D'^ c{x,y) is the identity matrix, and not merely constant. 

Our first contribution is the following theorem. If the cost function satisfies (B3), then 
the level sets of the c-convex potential appear convex from y, as was discovered independently 
from us by Liu |40) . and exploited by Liu with Trudinger and Wang |41) . Moreover, for a 
non-negatively cross-curved cost (B4), it shows that any c-convex potential appears convex 
from y £ V. Note that although the difference between the cost c{x, y) and the modified cost 
c{q, y) depends on y, they differ by a null Lagrangian c(x, y) which — being independent of 
y £ V — does not affect the question of which maps G attain the infimum (jl.ip . Having a 
function with convex level sets is a useful starting point, since it opens a possibility to apply 
the approach and techniques developed by Caffarelli, and refined by Gutieerrez, Forzani, 
Maldonado and others (see |28| |24|). to address the regularity of c-convex potentials. 

Theorem 4.3 (Modified c-convex functions appear level-set convex). Let c G C^(U x V) 
satisfy (B0)-(B3). If u = u^ ^ is c-convex on U, then u{q) = u{x{q))-\-c{x{q),y) has convex 
level sets, as a function of the cost exponential coordinates q £ Uy from y £ V . Moreover, 

u + Mc\q\ is convex, (4.5) 

where Mc := {(3~)'^\\Dl^c\\Loc(^u^v) + il3c)^\\Dxc\\L^{UxV)\\DlxDyc\\Lcx,(^Uy^V)- If in addition, 
c is non-negatively cross-curved (B4) then u is convex on Uy. In either case u is minimized 
c-i Qo if y £ d^u{x{qQ)) ■ Furthermore, u is c-convex with respect to the modified cost c{q,y) := 
c{x{q),y) — c{x{q),y) on Uy x V, and d^u{q) = d'^u{x{q)) for all q £ Uy. 

Proof. The final sentences of the theorem are elementary: c-convexity u = u^ ^ asserts 
u{x) = sup— c(x,y) — u^ (y) and u'^ (y) = sup —c{x{q),y) — u{x{q)) = u'^ (y) 

yeV q&Uy 



from ()3.2p . hence 



u{q) = sup -c{x{q),y) + c{x{q),y) -u" (y) 
ydV 

= sup -c{q,y) -u'^'iy), 

ydV 
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and d^u{q) = d'^u{x{q)) since all three suprema above are attained at the same y £ V. Taking 
y = y reduces the inequality u{q) + if {y) + c(g, y) > to u{q) > —vP (y) , with equality 
precisely if y € d^u{q). It remains to address the convexity claims. 

Since the supremum u{q) of a family of convex functions is again convex, it suffices to 
establish the convexity oi q £ Uy i — > —c{q,y) for each y GV under hypothesis (B4). For a 
similar reason, it suffices to establish the level-set convexity of the same family of functions 
under hypothesis (B3). 

First assume (B3). Since 

Dyc{q,y) = Dyc{x{q),y) := -q (4.6) 

we see that c-segments in Uy with respect to y coincide with ordinary line segments. Let 
q{s) = {l — s)qQ + sqi be any line segment in the convex set Uy. Define /(s, y) := —c{q{s), y) = 
—c{x{q{s)),y) + c{x(q(s)),y). Loeper's maximum principle (Theorem 13.31 above, see also 
Remark 13. 6p asserts f{s,y) < max{/(0,y), /(l,y)}, which implies convexity of each set 
{q G Uy I —c{q,y) < const}. Under hypothesis (B4), Theorem 13.51 goes on to assert 
convexity of s G [0, 1] i — > /(s, y) as desired. 

Finally, (14. 5 p follows from the fact that u is a supremum of cost functions and that 
\Dqqc\ < {f3^)'^\Dl^c\ + {/3~)^\Dxc\\Dl^Dyc\ (observe that Dqx{q) = —Dlyc{x{q),y)~^). D 

Remark 4.4. The above level-set convexity for the modified c-convex functions requires 
(B3) condition, as one can easily derive using Loeper's counterexample |42| . 



The effect of this change of gauge on the c-Monge- Ampere equation (13. Sp is summarized 
in a corollary: 

Corollary 4.5 (Transformed c-Monge- Ampere inequalities). Using the hypotheses and no- 
tation of Theorem \4-3[ if \d'^u\ G [A, A] C [0,cxd] on U' C U, then \d^u\ G [A/7^,A7~] on 
U'~ = -Dyc{U',y), where 7± = 7±(C/' x V) and pf = l3f{U' x V) are defined in (IO|l - (IOp . 
Furthermore, -^f := -if{U'y xV)< 7+7" and pf := (3f{UL x V) < /3+/3-. 

Proof. From the Jacobian bounds \det Dxq{x)\ G [^/^c^lc^] '^^ ^'1 '^^ ^^^ ■^"'i-^)/^c — 
^'^{q{X)) < 7+=Sf"(X) for each X C U'. On the other hand, Theoremgjasserts (9^u(g(X)) = 
d^u{X), so the claim \d'^u\ G [A/7^,A7~] follows from the hypothesis \d'^u\ G [A, A], by def- 
inition ()3.6p and the fact that q : U — > Uy from (14. ip is a diffeomorphism; see (Bl). The 
bounds 7g < ^^Ic ^'^^ Pc — Pt l^c follow from D'^yc{q,y) = D'^yC{x{q),y)Dqx{q) and 
Dqx{q) = -Dlyc{x{q),yy\ D 

Afflne renormalization 

We record here an observation that the c-Monge- Ampere measure is invariant under an affine 
renormalization. This is potentially useful in applications, though we do not use this fact for 
the results of the present paper. For an affine transformation L : R" — ?■ R", define 

n*(g) = |detL|-2/"u(Lg). (4.7) 
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Here det L denotes the Jacobian determinant of L, i.e. the determinant of the linear part of 
L. 

Lemma 4.6 (Affine invariance of c-Monge- Ampere measure). Assuming (BO)-(Bl), given 
a c- convex function u : Uy i — > R and affine bijection L : R" i — > R", define the renornialized 
potential u* by (j4.7p and renornialized cost 

c,iq, y) = I det L|-2/"c(L(/, L*y) (4.8) 

using the adjoint L* to the linear part of L. Then, for any Borel set Q C Uy, 

\du*\{L~^Q) = \ det L\-^\du\{Q), (4.9) 

\d^'u*\{L~^Q) = \detL\-^\d^u\{Q). (4.10) 



Proof From (f3l4|) we see p e du{q) if and only if | detL|~^/"L*p e du*{L~'^q), thus (|49]l 
follows from du*{L-'^Q) = | det L|-2/"L*(an(Q)). Similarly, since ([32]) yields {u*y*{y) = 
I det L|~^/"u'^*(L*y), we see y £ d^u{q) is equivalent to | detL|~^/"L*y G d'^'u*{L~^q) from 
(123D (and TheoremiS]), whence (9^*n*(L~iQ) = | det L|-2/"L*(0'^i(*(g)) to establish (ginD- 

D 

As a corollary to this lemma, we recover the affine invariance not only of the Monge- 
Ampere equation satisfied by u{q) — but also of the c-Monge- Ampere equation it satisfies 
— under coordinate changes on V (which induce linear transformations L on T^V and L* 
on TyV): for q £ Uy, 

d|9n*L,_i , d\du\, , , d\d^*u*\^^, , dld^ul , , 

^^(^- «) = ^^(«) and ^^^^(^- «) = ^^('?)- 

5 Strongly c-convex interiors and boundaries not mixed by d^u 

The subsequent sections of this paper are largely devoted to ruling out exposed points in Uy 
of the set where the c-convex potential from Theorem 14.31 takes its minimum. This current 
section rules out exposed points on the boundary of Uy. We do this by proving an important 
topological property of the (multi-valued) mapping d'^u C U xV. Namely, we show that the 
subdifferential d'^u maps interior points of spt |(9'^m| C U only to interior points of V, under 
hypothesis (j3.8p , and conversely that d'-^u maps boundary points of U only to boundary points 
of V. This theorem may be of independent interest, and was required by Figalli and Loeper 
to conclude their continuity result concerning maps between two dimensional domains which 
optimize (B3) costs |20) . 

This section does not use the (B3) assumption on the cost function c G C^{U x V), but 
relies crucially on the strong c-convexity (B2)s of its domains U and V (but importantly, not 
of spt l^^^ul). No analog for Theorem 15.11 was needed by Caffarelli to establish C^'" regularity 
of convex potentials u{x) whose gradients optimize the classical cost c{x,y) = —{x,y) [8], 
since in that case he was able to take advantage of the fact that the cost function is smooth 
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on the whole of R" to chase potentiahy singular behaviour to infinity. (One general approach 
to showing regularity of solutions for degenerate elliptic partial differential equations is to 
exploit the threshold-hyperbolic nature of the solution to try to follow either its singularities 
or its degeneracies to the boundary, where they can hopefully be shown to be in contradiction 
with boundary conditions; the degenerate nature of the ellipticity precludes the possibility of 
purely local regularizing effects.) 

Theorem 5.1 (Strongly c-convex interiors and boundaries not mixed by d'^u). Let c satisfy 
(BO)-(Bl) and u = u'^ '^ he a c-convex function (which implies d'^u{U) = V), and A > 0. 

(a) If Id'^ul > X on X C U and V is strongly c* -convex with respect to X, then interior 
points of X cannot he mapped hy d'^u to houndary points ofV: i.e. {X x dV) n d^u C 
{dX xdV). 

(h) If \d^u\ < A on U, and U is strongly c-convex with respect to V , then houndary points 
of U cannot he mapped by d^u into interior points ofV: i.e. dU x V is disjoint from 
d'^u. 

Proof. Note that when X is open the conclusion of (a) implies d'^u is disjoint from X x dV . 
We therefore remark that it suffices to prove (a), since (b) follows from (a) exchanging the 
role X and y and observing that \d'^u\ < A implies \d'^ u^ \>1/ K as in Lemma l3.ir d). 

Let us prove (a). Fix any point x in the interior of X, and y E d^u{x). Assume by 
contradiction that y G dV. The idea of proof is summarized in Figure [1] We first fix 
appropriate coordinates. At {x,y) we use (BO)-(Bl) to define cost-exponential coordinates 
{p,q) I — > {x{q),y{p)) by 

p= -D,c{x,y{p)) + D^c{x,y) G r|(^) 

q= DlyC{x,yr\Dyc{x{q),y)-Dyc{x,y)) eT.j,{U) 

and define a modified cost and potential by subtracting null Lagrangian terms: 

c{q,p) := c{x{q),y{p)) - c{x{p),y) - c{x,y{p)) 
u{q) := u{x{q))+c{x{q),y). 

Similarly to Corollary 14.51 \d'^u\ > A := A/(7^7~), where 7^ denote the Jacobian bounds 
(j4.4p for the coordinate change. Note {x,y) = (x(0),y(0)) corresponds to {p,q) = (0,0). 
Since c-segments with respect to y correspond to line segments in Uy := —Dyc{U,y) we see 
Dpc{q, 0) depends linearly on q, whence D^ c{q, 0) = 0; similarly c*-segments with respect 
to X become line segments in the p variables, Dqc{0,p) depends linearly on p, D^ c{0,p) = 0, 
and the extra factor D"^ c{x,y)~^ in our definition of x{q) makes —Dpc{0,0) the identity 
matrix (whence q = —Dpc{0,q) and p = —Dqc{p,0) for all q m Uy = x~^{U) and p in 
Vx := y~^{V)). We denote Xy := x~^{X) and choose orthogonal coordinates on U which 
make — e^ the outer unit normal to Vj. C T|C/ at p = 0. Note that Vx is strongly convex by 
hypothesis (a). 



15 





Figure 1: If d'^ii sends an interior point onto a boundary point, by c-monotonicity of d'^u the small 
cone Eg^i; (with height e and opening 9) has to be sent onto Eg ^^^ n Vj. Since ^""(-Bg.e) ^ 0"^^ for 
e > small but fixed, while ^^^{E'g ^^^ n 14) < 0"+^ (by the strong convexity of Vx), we get a 
contradiction as — )■ 0. 



In these variables, as in Figure (TJ consider a small cone of height e and angle around 
the — e„ axis: 



-E'e.e : = 



GR" 



< sin0, |g| < e 



Observe that, if 0,e are small enough, then Eg^^ C Xy, and its measure is of order e"0" ^. 
Consider now a slight enlargement 



E', 



Cos 



P 



{P,Pn)€'R.^\Pn<0\p\+Coe\p\^}, 



of the polar dual cone, where e will be chosen sufficiently small depending on the large 
parameter Cq forced on us later. 

The strong convexity ensures Vx is contained in a ball Bpi{Ren) of some radius R > 1 
contained in the half-space Pn ^ with boundary sphere passing through the origin. As long 
as CqE < {6R)~^ we claim E'q ^j ^ intersects this ball — a fortiori Vx — in a set whose volume 
tends to zero like 0"+^ as — )■ 0. Indeed, from the inequality 

Pn<eyJ\P\^+pl + \\P\ + \pn 

b 6 

satisfied by any (P,p„) G E'q(.^^ n BniRin) we deduce pi < \P\'^{1 + 96'2)/(2 - 96^2), i.e. 
Pn < \P\ li 9 is small enough. Combined with the further inequalities 



2R 



<Pn< OV\P\^+pI + Coe\P\^ + Coepl 



(the first inequality follows by the strong convexity of Vx , and the second from the definition 
of E'oc^), this yields |P| < 60\/2 and pn < 0{9^) as 6 ^ 0. Thus ^""{E'g^^,^ n V^) < C0"+i 
for a dimension dependent constant C, provided Cq£ < {6R)~^ . 
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The contradiction now will come from the fact that, thanks to the c-cyclical monotonicity 
of d'^u, if we first choose Cq big and then we take e sufficiently small, the image of all q G Eg^^ 

by d'^u has to be contained in E'q (j ^ for 9 small enough. Since d^uiXy J C Vj this will imply 

which gives a contradiction as — )■ 0, for e > small but fixed. 

Thus all we need to prove is that, if Cq is big enough, then d^u{Eg^£) C Eg q ^ for any e 
sufficiently small. Let q G Eg^^ and p £ d^u{q). Notice that 

ds dtDgpC{sq,tp)[q,p] 
Jo Jo 

= c{q,p)+c{0,0)-c{q,0)-c{0,p) 

< (5.1) 

where the last inequality is a consequence of c-monotonicity of d'^u, see for instance |54| 
Definitions 5.1 and 5.7]. Also note that 

DgpC{sq,tp) =Dg c{0,tp) + / ds'D^gpC{s'q,tp)[q] 

Jo 

=Z?2c(o,o)+ / dt'DlppC{0,t'p)[p] 
Jo 

ps ps pt 

+ / ds'D^ggpC{s'q,0)[q]+ / ds' / dt'D^gppC{s'q,t'p)[q,p] 
Jo Jo Jo 

= -In+ I ds' [ dt'DlgppC{s'q,t'p)[q,p] (5.2) 

Jo Jo 

since D^ c{0,t'p) and D^ c{s'q,0) vanish in our chosen coordinates, and —Dpc{0,0)= In 
is the identity matrix. Then, plugging (15. 2p into (15. ip yields 

pi pi ps pt 

-{q,p) <- / ds dt ds' / dt'D^gppC{s'q,t'p)[q,q,p,p] 
Jo Jo Jo Jo 

< Co\q\'\p\ 



for constant Cq = su]ptgp\^x-xVs: W^qgpp'^W- Since the term —D^^ppC is exactly the cross- 
curvature (j2.ip . its tensorial nature implies Cq depends only on ||c||(^4(f/xy) and the bi- 
Lipschitz constants /3^ from (|4.3p . 

From the above inequality and the definition of Eg^^ we deduce 

Pn = {p, en + rr) - (p^ rr) < ^1^1 + ^o^IpP 
\q\ \q\ 

so p £ E'g pp^ as desired. D 
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6 Alexandrov type estimates for c- convex functions 

In this section we prove the key estimates for c-convex potential functions u which will 
eventually lead to the Holder continuity and injectivity of optimal maps. Namely, we extend 
Alexandrov type estimates commonly used in the analysis of convex solutions to the Monge- 
Ampere equation with c{x,y) = —{x,y), to c-convex solutions of the c-Monge- Ampere for 
general (B3) cost functions. These estimates, Theorems 16.21 and 16.111 are of independent 
interest: they concern sections of u, i.e. the convex sub-level sets of the modified potential 
u of Theorem 14. 3| and compare the range of values and boundary behaviour of u on each 
section with the volume of the section. We describe them briefly before we begin the details. 
Fix p G R", pq G R, and a positive symmetric matrix P > 0. It is elementary to see that 
the range of values taken by the parabola u{q) = q^Pq+p-q+po on any non-empty sub-level 
set Q = {u < 0} is determined by det P and the volume of Q: 

I min u(q)r = ( 5"/'^l ^ det P. (6.1) 

Moreover, the parabola tends to zero linearly as the boundary of Q is approached. Should 
the parabola be replaced by a convex function satisfying A < det D^u < A throughout a fixed 
fraction of Q, two of the cornerstones of Caffarelli's regularity theory are that the identity 
(16. ip remains true — up to a factor controlled by A, A and the fraction of Q on which these 
bounds hold — and moreover that u tends to zero at a rate no slower than distgQ{q)^'"' as 
q — )• dQ. The present section is devoted to showing that under (B3), similar estimates hold 
for level-set convex solutions u of the c-Monge Ampere equation, on sufficiently small sub- 
level sets. Although the rate distgQ(g)^'^ we obtain for decay of u near dQ is probably 
far from optimal, it turns out to be sufficient for our present purpose. 

We begin with a Lipschitz estimate on the cost function c, which turns out to be useful. 

Lemma 6.1 (Modified cost gradient direction is Lipschitz). Assume (B0)-(B2). Fixy G V. 
For c G C^\Uy X V) from Definition \4-l\ and each q,q G Uy and fixed target y (^ V, 

I - Dgc{q,y) +Dgc{q,y)\ < —\q-q\ \Dgc{q,y)\, (6.2) 

where Sc is given by e^^ = '^iPt)'^il^c)^\\^xxy'^\\L°°{UxV) ^'^ ^^^ notation (j4.3p . If y ^ y so 
neither gradient vanishes, then 



Dqc{q,y) Dqc{q,y) 



\Dqc{q,y)\ \Dqc{q,y)\ 



<-\q-q\. (6.3) 



Proof. For fixed q G Uy introduce the c-exponential coordinates p{y) = —Dqc{q,y). The 
bi-Lipschitz constants (|4.3p of this coordinate change are estimated by f]~ < /3^/3~ as in 
Corollary 14.51 Thus 



dist(y, y) <l3-\- Dqc{q, y) + Dqc{q, ; 



<PtPc\Dqc{q,y)V 
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where c{q,y) = from Definition 14.11 has been used. Similarly, noting the convexity (B2) of 



V, 



p{V), 



Dqc{q,y) + Dqc{q,y)\ 



Dqc{q, y) + Dqc{q, y) + Dqc{q, y) - Dqc{q, ', 



< \\DtD„c 



■qqDpc\\^^^^_^y_)\q - q\\p{y) - p{y)\ 



< \\DqqDyc\\L^^u.^v)il3c f^c) |g-9|dist(y,y) 

Then dEl follows since \D^^gDyc\ < jC/^c^)' + /3+(/5c^)'))I^L^.c| < 2/3+(/3,-)3|I)2^Z?,c| . 
(The last inequality follows from P^/3^ ^1-) 

Finally, to prove (|6.3p we simply use the triangle inequality and ()6.2p to get 



Dqc{q,y) Dqc{q,y) 



\Dqc{q,y)\ \Dqc{q, 



< 



Dqc{q,y) Dgc{q,y) 



\Dqc{q,y)\ \Dqc{q,y)\ 
\Dqc{q,y)\ + \Dqc{q,y)\\ 



\Dqc{q,y)\ 



+ 



+ 



Dqc{q,y) Dyc{q,y) 



\Dqc{q,y)\ \Dqc{q,y)\ 

Dqc{q,y) +Dqc{q,y)\ 
\Dqc{q,y)\ 



2 , 

< — 



D 



6.1 Alexandrov type lower bounds 

In this subsection we prove one of the key estimates of this paper, namely, a bound on the 
infimum on a c-convex function inside a section in terms of the measure of the section and 
the mass of the c-Monge- Ampere measure inside that section. The corresponding result for 
the affine cost function is very easy, but in our case the unavoidable nonlinearity of the cost 
function requires new ideas. 

Theorem 6.2 (Alexandrov type lower bounds). Assume (B0)-(B3), define c G C^{Uy x V) 
as in Definition \4-l\ and let -)~ = 7^ (Q x V) he as in ( f^.^p . Let u : Uy i — > H be a c-convex 



IS convex. 



function as in Theorem \4-3[ and let Q := {u < 0} C Uy. Note that Q C R" 
Let Ec > small be given by Lemma \6.1\ Finally, let E be the ellipsoid given by John's 
Lemma (see (j3.9p ). and assume there exists a small ellipsoid Eg = xq + 5E C Q, with 
6 < min{l,ec/(4diam(i?))} and {d'^ul > A > inside Eg. We also assume that Eg <Z B C 
4B C Uy, where B, 4B are a ball and its dilation. Then there exists a constant C{n), 
depending only on the dimension, such that 



^'\Qf < Cin] 



7.- 



52" A 



infiir. 



(6.4) 



The proof of the theorem above, which is given in the last part of this subsection, relies 
on the following result, which will also play a key role in Section [9] (see (j9.4p ) to show the 
engulfing property and obtain Holder continuity of optimal maps. 
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Figure 2: The quantity ^^-rr^ represents tlie distance between the origin and the supporting hyper- 
plane orthogonal to v. 



Lemma 6.3 (Dual norm estimates). With the same notation and assumptions as in Theorem 
\6.^ let J(f d Q he an open convex set such that diam(^) < ec/4. We also assume that there 
exists a hall B such that J^ C B C 4:B C Uy, where AB denotes the dilation of B by a factor 
4 with respect to its center. (This assumption is to locate all relevant points inside the domain 
where the assumptions on c hold.) Up to a translation, assume that the ellipsoid associated 
to .J€ by John's Lemma fsee \3.9\) is centered at the origin 0, and for any p G (0, 1) let pj^ 
denote the dilation of J€ with respect to the origin. Moreover, let \\ ■ \\*-^ denote the "dual 
norm" associated to J(f , that is 

\\v\\*j(f := sup w-v. (6-5) 

■u;G,X 

Then, for any p < 1 setting C^{n, p) := /i_\i implies 

\\-Dgc{q,y)\\*^<C^{n,p)\miu\ y q e p,^ , y e d^u{p,^). (6.6) 

Before proving the above lemma, let us explain the geometric intuition behind the result: 
for V G R" \ {0}, let H^ be the supporting hyperplane to J(f orthogonal to v and contained 
inside the half-space {g | g • f > 0}. Then 

dist(0, Hy) = sup w ■ -r—r = — r^, 

w<E.T \V\ \V\ 

see Figure [21 and Lemma 16.31 states that, for all v = —Dqc{q,y), with q € pj^ and y G 
d^u{pJ(r), 

dist{0, Hy)\v\ < C*(n, p) I inf-ul. 

Hence, roughly speaking, (j6.6p is just telling us that the size of gradient at a point inside p^ 
times the width of ^ in the direction orthogonal to the gradient is controlled, up to a factor 
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Figure 3: It is easily seen that \qo — q\ controls dist(0, H^,), and using that the gradient of niy does 
not vary much we deduce that the latter is close to dist(0, H^) in terms of |go — <z|- 



Ct:{n,p), by the infimum of u inside ^ (all this provided the diameter of ^ is sufficiently 
small). This would be a standard estimate if we were working with convex functions and 
c was the standard quadratic cost in R", but in our situation the proof is both subtle and 
involved. 

To prove (16. 6p let us start observing that, since \w\ < diaiQ{J^) for w G J(f, the following 
useful inequality holds: 

II^ILV < diam(jr) |f I Vu G R". (6.7) 

We will need two preliminary results. 

Lemma 6.4. With the same notation and assumptions as in Lemma \6.3[ let q G p-^ , y G V, 
and let niy be a function of the form niy := —c{-,y) + Cy for some constant Cy G R. Set 
V := —Dqiriy^q), assume that J(f C {niy < 0}, and let qo denote the intersection of the 
half-line i1 '■= q + R~^v = {q + tv\t > 0} with {niy = 0} (assuming it exists in Uy). Define 



qo 



Q + 



diam(iJ^) 



*/ ko — ^l < diamJ^, 
V if ko — ?! ^ diam(^) or go does not exist in Uy, 



(Notice that, by the assumption ^ d B C 4i? <Z Uy, we have qo G Uy.) Then 



2\qo-q\>{l-p)- 



\X 



(6.8) 



Proof. In the case % = q -\- diaui{J^)v , the inequality (j6.8p follows from (|6.7p . Thus, we can 
assume that qo exists in Uy and that % = go- Let us recall that {niy < 0} is convex (see 
Theorem 14. 3p . Now, let H^ denote the hyperplane tangent to {niy = 0} at qo. Let us observe 
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that Hi, is orthogonal to the vector vq := —Dqmy{qo), see Figure [3l Since J^ C {rriy < 0}, 
we have Hy D J(f = 0, which imphes 

dist(0, Hyj > sup w ■ - — j- = — — f-. 

wex \vo\ \vo\ 

Moreover, since q E p^ and H^ J^ = 0, we also have 

dist{q,Hy) > (1 -p)dist(0,^„). 
Hence, observing that go G H^, we obtain 



|go - q\> dist(g, Hy) > (1 - p) 



\mx 



Now, to conclude the proof, we observe that (j6.3p applied with vq = —Dqc{qo,y) and f 
—Dqc{q,y), together with (16. 7p and the assumption diam(J^) < ec/4, implies 



^^0 

l^ol 



< diam(^) 



X 






2diam(=J^) I go 
< ko - g| < ; 



Combining all together and using the triangle inequality for || • ||^ (which is a consequence 
of the convexity of J(f) we finally obtain 



V 

\v\ 


* 

< 

X 


\vo\ 


* 

+ 

X 






X 



< 



< 



m 



i-p 
2^0 - q 



+ 



1 



as desired. 



D 



Lemma 6.5. With the same notation and assumptions as in Lemmata \6.3\ and \6.4\ fix 

I I 

q' G P'^ , o,nd let 1% denote the half line H.% := q' + R+w = {q' + tv \t > 0}. Denote hy q^ 
the intersection of 1% with {my = 0} (assuming it exists in Uy). Then, 



> 



2n 



m - q\ 



Proof Let ^g^/ and ig^gi denote the lines passing through q,q' and goi^o respectively, and 
denote by qi := igg' H i^^^g' their intersection point (see Figure [4]); since the first four points 
lie in the same plane, we can (after slightly perturbing q or q' if necessary) assume this 
intersection exists and is unique. Note that in the following we will only use convexity of J^, 
and in particular we do not require qi € Uy. 

Let us now distinguish two cases, depending whether |g— gi| < |g' — gi| or |g— gi| > |g' — gi|- 
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Figure 4: In this figure we assume go — Qo- Using similarity, we only need to bound from below the 



ratio 



— and the latter is greater or equal than '? _ ? , . An easy geometric argument allows to 



I9-9 
control this last quantity. 



I9-91I 



• k ~ 91! ^ W ~ Qi\- III this case we simply observe that, since ^0 
parallel, by similarity 

W -qi\ _ Wo-q'\ 



q and 



fn - q' are 



1 < 



\q-qi\ 



m 



and so the result is proved since 1 > — - 



2n 



• k ~ ^il ^ \q ~ Qi\'- III this second case, we first claim that qi ^ J^ . Indeed, we observe 
that since qo ~ Q and gg — q' are parallel, the point qi cannot belong to the segment joining 
qo and q'^. Hence, since ^0 £ {"t-j/ < 0}, (/q € {niy = 0}, and Jf C {my < 0} by assumption, 
by the convexity of the set {niy < 0} we necessarily have qi outside {niy < 0}, thus outside 
of J^, which proves the claim. 

Thus, we can find the point qi obtained by intersecting d,J^ with the segment going from 
q' to qi, and by the elementary inequality 



a + c a 
b + c ~ 1) 



VO<a< 6, c> 0, 



we get 



Qil 



> 



Qi\ 



\Q-Qi\ iQ-qil 
Now, to estimate the right hand side from below, let L be the affine transformation provided 
by John's Lemma such that Bi C L~^(^) C Bn- Since the points q,q',qi are aligned, we 
have 



Qil 



\L-^q' 



L^'qil 



L-^i\ 



\q-qi\ \L ig 
Now, since L~^q,L~'^qi G Bn, we immediately get 

\L-^q- L~^qi\ < 2n. 
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On the other hand, since L ^q' G p(^L ^{J^)) while L ^qi G d[L ^(^)), it follows from 
Bi C L-i(^) that 

\L-\' -L-'^qil > 1-/3. 

To see this last inequality, for each point w £ 5(L~^(^)") consider the convex hull Cw C 

L^^(^) of {w} U Bi. Then, centered at the point pw £ dip(^L^^(Jif)) I there exists a ball 

of radius 1 — p contained in C^, thus in L~^{J^). This shows that the distance from any 
point in p(L^^(^)) to (9(L~^(^)) is at least 1 — p, proving the inequality. 
From these inequalities and using again similarity, we get 

Wo - q'\ ^ |g^-gil ^ 1-p 
\Qo-q\ \q-qi\ ~ 2n 

concluding the proof. D 

Proof of Lemma 16.31 Fix v = —Dqc{q,y), with q G pJ^, q' G p^ and y G d'^u{q'), and 
define the function my{-) := —c{-,y) + c{q',y) + u{q'). Notice that since y G d'^u{q'), the 
inclusion Q = {n < 0} C {niy < 0} holds. 

Recall the point qQ from Lemma 16.51 (when it exists in Uy), and the point qo G Uy from 
Lemma [6. 4[ Define, 



„, I q'o if ko ~ ^'1 — diam^. 



'^ ' I 5' H — '^"f I ' f if I^Q — (/'I > diam(J??') or q'^ does not exist in Uy, 



Notice that by the assumption ,J^ C B C AB C Uy, the points qQ, cj'q belong to Uy. Let 
[<?')9c)] '^ ^y denote the segment going from q' to ^q. Then, since ^q G {my < 0} and u is 
negative inside J^ C Q, we get 

I inf m| > —u{q') > —u{q') + n(^o) 

>-c{q'o,y) + ciq',y)= [ {-D,c{w,y),^) dw (6.9) 

Since q, g' G =>^ and \qQ — q'\ < diam(^), we have 

!(? — w| < 2diam(J^) < ec/2 Vw G [(7', (^q]- 
Hence, recalling that —Dqc{q,y) = v and applying Lemma l6.H we obtain 

{-Dqc{w,y),—)>— Vu;G[g',go]. 
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Applying this and Lemmas 16.41 and 16.51 to (j6.9p , we finally conclude 
I inf n| > min {diam(J^), |gg — q'\} — (from (|6.9p ) 

> min < diam(J^), I^q ~ ^I f tt (from Lemma |6. 5 



2 

f { '\ \ 2 1 1 1 1 * "^ 

> min < diam(J^), —^ > ^—^ (from Lemma [67 

[ 4n |t;| J '^ 



2 



• \^ A- ( ^w I iLli^ii II* 
= mm<^ -diam(^)|u|, — — — \\v\\j^ 

>mm|-||7;||,^, — ^;^Flljr| (from ([6 

-(1-/^)^1,1. 
8n " "-^ 

as desired, completing the proof. D 

Proof of Theorem [Ol Set h := \ infg u\. Since ^'"(Q) < n"^"(£;), it is enough to show 

X^^iEf < C(n)7,-^. (6.10) 

So, the rest of the proof is devoted to (|6.10p . Let Eg = xq + 6E be defined as in the statement 
of the proposition. With no loss of generality, up to a change of coordinates we can assume 
that Es is of the form {q \ ^^ a^qt^ < !}• Define 

Co = {qo G TqQ I (/o = -Dgc{0,y), y G d'u{lEs)} . 

Notice that 

^^{d'u{lEs))<^^^^{Co) (6.11) 

(see (14.41)). Now, let us define the norm 



1911b, 



\ 



n 



-\= sup V • g. 



i=l « 

Since J(f = Eg satisfies the assumptions for Lemma [6.31 choosing p = 1/2 we get 

\\qo\\*E, < ^"inh for all q^ G Cq. (6.12) 

Let (^h : R" — )• R" denote the linear map ^h{x) := 32n/i (aixi, . . . , anXn) ■ Then from (|6.12p 
it follows 

Co C ^h{Bi). (6.13) 
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Since 



^-{^k{B,))) < dn)^^^, (6.14) 



combining (j6.1ip . (j6.13p . and (j6.14p . we get 

h" 



^^{d'u{lEs))<C{n)j^ 



^""{Es 



As ^"(9S(i^5)) > A2-"if"(E5) and ^"(^5) = <5"^"(^), we obtain the desired conclu- 
sion. D 

6.2 Bounds for c-cones over convex sets 

We now progress toward the Alexandrov type upper bounds in Theorem 16.111 In this sub- 
section we construct and study the c-cone associated to the section of a c-convex function. 
This c-cone — whose entire c-Monge- Ampere mass concentrates at a single prescribed point 
— plays an essential role in our proof of Lemma 16.101 

Definition 6.6 (c-cone). Assume (B0)-(B3), and let u : Uy 1 — > R be the c-convex function 
with convex level sets from Theorem \4-3\ Let Q denote the section {u < 0}, fix q £ intQ, 
and assume u = on dQ and Q C Uy. The c-cone h^ : Uy 1 — > R generated by q and Q with 
height —u{q) > is given by 

h^iq) := sup{-c(g, y) + c{q, y) + u{q) \ -c{q, y) + c{q, y) + u{q) < on dQ}. (6.15) 

y&V 

Notice the c-cone h'^ depends only on the convex set Q C Uy, q G intQ, and the value 
■u(g), but is otherwise independent of u. Recalling that c{q,y) = on Uy, we record several 
key properties of the c-cone: 

Lemma 6.7 (Basic properties of c-cones). Adopting the notation and hypotheses of Definition 
\6.6^ let h^ : Uq 1 — > R be the c-cone generated by q and Q with height —u{q)> 0. Then 

(a) h'^ has convex level sets; furthermore, it is a convex function if (B4) holds; 

(b) h%q) > h%q) = u{q) for all q G Q; 

(c) h^ = on dQ; 

(d) d^h^iq) cd^u{Q). 

Proof. Property (a) is a consequence of the level-set convexity of q 1 — > —c{q, y) proved in 
Theorem 14. 3| or its convexity assuming (B4). Moreover, since —c{q,y)-\-c{q,y)-\-u{q) = u{q) 
for all q £ Uy, (b) follows. 

For each pair qo S dQ and yo S d'^u{qQ), consider the supporting mountain mQ{q) = 
—c{q,yo) + c{qo,yo), i.e. mo(go) = = u{qo) and tuq < u. Consider the c-segment y{t) 
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connecting y(0) = yo and y(l) = y inV with respect to go- Since —c{q,y) = 0, by continuity 
there exists some t G [0, 1[ for which mt{q) := —c{q,y{t)) + c{qo,y{t)) satisfies mt{q) = u{q). 
From Loeper's maximum principle (Theorem 13.31 above), we have 

mt < max[mo,-c(-,y)+c(go,y)] = max[mo,0], 

and therefore, from itiq < u, 

mt < on Q. 

By the construction, mt is of the form 

-ci-,yit))+c{q,y{t))+u{q), 

and vanishes at go- This proves (c). Finally (d) follows from (c) and the fact that h'^{q) = u{q). 
Indeed, it suffices to move down the supporting mountain of h'^ at q until the last moment 
at which it touches the graph of u inside Q. The conclusion then follows from Loeper's local 
to global principle, Corollarv 13.41 above. D 

The following estimate shows that the Monge- Ampere measure, and the relative location 
of the vertex within the section which generates it, control the height of any well-localized 
c-cone. Together with Lemma iG.Tf d). this proposition plays a key role in the proof of our 
Alexandrov type estimate (Lemma 16. lOh . 

Proposition 6.8 (Lower bound on the Monge- Ampere measure of a small c-cone). Assume 
(B0)-(B3), and define c G C^^UyXV^ as in Definition \4-l\ Let Q C Uy be a compact convex 
set, and h^ the c-cone generated by q £ mtQ of height —h'^{q) > over Q. Let 11"'", H" be two 
parallel hyperplanes contained in T^V\ Q and touching dQ from two opposite sides. We also 
assume that there exists a ball B such that Q G B C AB C Uy. Then there exist ffc > small, 
depending only on the cost (and given by Lemma \6.1\) . and a constant C{n) > depending 
only on the dimension, such that the following holds: 
//diam((5) < e^ ■.=£c/C{n), then 

mqT < C(^) "^"^idist(g,n+),dist(g,n-)} ^^^,^^^,^^^„^^^^ ^g^^g^ 

where 1-^+ denotes the maximal length among all the segments obtained by intersecting Q with 
a line orthogonal to H"'". 

Proof. We fix q£ int Q. Let H*, i = 1, - - - n, (with H^ equal either H^ or H^) be hyperplanes 
contained in Ti^V \ Q ~ R" \ Q, touching dQ, and such that {n+,n^, . . . ,n"} are all 
mutually orthogonal (so that {11", H^, . . . ,n"} are also mutually orthogonal). Moreover we 
choose {n^, . . . ,n"} in such a way that, if ■7r^(Q) denotes the projection of Q on H^ and 
J^"'~^ {tt^ (Q)) denotes its (n — l)-dimensional Hausdorff measure, then 

n 

C7(n) jr"-i(^i(Q)) > Yl dist(g, ff), (6.17) 



i=2 
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for some universal constant C{n). Indeed, as vr^(Q) is convex, by Lemma [3.71 we can find an 
ellipsoid E such that E C vr^(Q) C (n — 1)E, and for instance we can choose {H^, . . . ,11'"} 
among the hyperplanes orthogonal to the axes of the ellipsoid (for each axis we have two 
possible hyperplanes, and we can choose either of them). 

Each hyperplane 11* touches Q from outside, say at q^ G T^V. Let pi € TyV be the 
outward (from Q) unit vector at g* orthogonal to W. Then SiPi G dh'^{q^) for some Sj > 0, 
and by Corollary 13.41 there exists yi G d^hP{q^) such that 

-Dqc{q\yi) = SiPi. 

Define yi{t) as 

-Dqc{q\yi{t)) = tSiPi, 

i.e. yi{t) is the c-segment from y to yi with respect to q^. As in the proof of Lemma 16.7( c). 
there exists < ij < 1 such that the function 

^y^{t^)i■) '■= -c{-,yi{ti)) + c{q,yi{ti)) + h'^iq) 

satisfies 

''^t/iiti) !^ on Q with equality at g*, (6.18) 

see Figure[5l By the definition of h'^, (I6.18J) implies yi{ti) £ d'^h^{q) Pi d'^h'^{q^), 

-Dqc{q,yi{ti)) G dh^'iq) and tiSiPi = -Dgc{q\yi{ti)) G dh'^iq'). 

Note that the sub-level set Sy.u.\ := {z G f/y | T^yJti) ^ 0} is convex. Draw from q a half-line 
orthogonal to IT*. Let g* be point where this line meets with the boundary dSy.u.y By the 
assumption Q C B C AB C Ug, we see q^ G Uy. By convexity of Sy.u.\, we have 

|g-gi < dist(g,rf). (6.19) 

Let diamQ < 6n£c for some small constant < (5„ < 1/3 to be fixed. By (|6.2p and the trivial 
inequality dist(g, H*) < diamQ, we have 

\-Dqc{q,y,{ti))+Dqc{iq+{l-^)q',yi{ti))\ < -\q - q'\\ - Dqc{q,yi{U)\ 

Ec (6.20) 

< 5n\ - DqC{q,yi{ti))\ 

for all < £ < 1. Therefore, we see 

\h^{Q)\ = \c{q,yt{U))-c{q\yi{ti))\ 

[ -D,c{£q + (1 - e)q\ y,{ti)) ■ {q - q')M 
'0 

<{l + 6n)\-Dqc{q,y^{U))\\q-q'\ (by ([620])) 

< (1 + <5„)| -Dqc(g,y,(t,))|dist(g,tf) by ^M)- 
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Figure 5: The dotted line represents tire graph of rriy^ := ~c{-,yi) + c(g, j/i) + h'^{q), while the dashed 
one represeirts the graph of 'my.(^t.^ := —c{-,yi{ti)) + c{q,yi{ti)) + h'^{q). The idea is that, whenever 
we have niy- a supporting function for h'^ at a point g* € dQ, we can let y vary continuously along 
the c-segment from y to yi with respect to g*, to obtain a supporting function rUy.uA which touches 
h'^ also at q as well as 5'. 



Thanks to this estimate it follows 



Dqc{q,yi{ti))\ > 



\h%q)\ 



2dist(g,m)' 
Moreover, similarly to (j6.20p . we have 

I - Dqc{q\yi{ti)) + Dgc{q,yi(ti))\ < 5„| - Dqc{q,yi{ti))\ 



(6.21) 



(6.22) 



Since the vectors {—Dqc{q^,yi{ti))}"^'^ are mutually orthogonal, (I6.22p implies that for 5n 
small enough the convex hull of {—Dqc{q,yi{ti))}^^^ C dh'^{q) has measure of order 

n 

ni -^gc(g,yi(ti))|. 

i=l 



Thus, by the lower bound ()6.2ip and the convexity of dh'^{q), we obtain 

^"(5/i'=(g)) > C7(n)=^— ^^L^. 
nr=idist(g,m) 

Since 11^ was either 11+ or H^, we have proved that 

n 

|/i^(Q^)|"<C7(n)|5/i^~|({g})min{dist(g,n+),dist(g,n-)}JJdist(g,rf). 

i=2 

To conclude the proof, we apply Lemma 16.91 below with Q' given by the segment obtained 
intersecting Q with a line orthogonal to 11'^. Combining that lemma with (|6.17p . we obtain 

n 

i=2 
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Figure 6: The volume of any convex set always controls the product (measure of one slice) • (measure 
of the projection orthogonal to the slice) . 



and last two inequalities prove the proposition (taking C{n) > l/5n larger if necessary). D 

Lemma 6.9 (Estimating a convex volume using one slice and an orthogonal projection). Let 
Q be a convex set in R" = R" x R" . Let tt', tt" denote the projections to the components 
R" , R" , respectively. Let Q' be a slice orthogonal to the second component, that is 

Q' = {7r")-\x") n Q for some x" G 7r"{Q). 

Then there exists a constant C{n), depending only on n = n' + n" , such that 

C7(n)^"(Q) > ^"'(Q')^""(^"(Q)), 
where .M"^ denotes the d-dimensional Hausdorff measure. 



Proof. Let L : R" — ?• R" be an affine map with determinant 1 given by Lemma 13.71 such 
that Br C L['k"{Q)) C Bnii^ for some r > 0. Then, if we extend L to the whole R" as 
L{x',x") = {x',Lx"), we have ^"(L(Q)) = ^"(Q), jr"'(L(Q')) = -^""'{Q'), and 

jr""(7r"(Z(Q))) = jr""(L(7r"(Q))) = ^""(7r"(Q))- 

Hence, we can assume from the beginning that Br C tt"{Q) C Bn^'r- Let us now consider 
the point x" , and we fix an orthonormal basis {ei, . . . ,e„//} in R" such that x" = cei for 
some c < 0. Since {rei, . . . ,re„//} C tt"{Q), there exist points {xi, . . . ,x„"} C Q such that 
'7r"{xi) = rei. Let C denote the convex hull of Q' with xi, and let V denote the {n' + 1)- 
dimensional strip obtained taking the convex hull of R" x {x"} with xi. Observe that 
C C V\ so 



^"'+i(C7') = — ^dist(xi,R"' X {x"})jr"'(Q') > — ^_Jr"'(Q')• 
n' + 1 n' + 1 



(6.23) 
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We now remark that, since 7r"(xj) = rej and e, _L V' for i = 2, . . . , n" , we have dist(xi, V') = r 
for all i = 2, . . . ,n". Moreover, if yi £ V' denotes the closest point to Xj, then the segments 
joining Xj to yi parallels e^, hence these segments are all mutually orthogonal, and they are 
all orthogonal to V' too. From this fact it is easy to see that, if we define the convex hull 

C := co{x2,...,Xn",C'), 

then, since \xi — yi\ = r for i = 2, . . . ,n", by (j6.23p and the inclusion tt"{Q) C Bn"r C R" 
we get 

n\ n\ 

This concludes the proof, as C C Q. D 

6.3 Alexandrov type upper bounds 

The next Alexandrov type lemma holds for localized sections Q of c-convex functions. 

Lemma 6.10 (Alexandrov type estimate and lower barrier). Assume (B0)-(B3), and let 

u : Uy I — 7- Yi he a c-convex function from Theorem \4-S\ Let Q denote the section {u < 0} C 
Uy, assume u = on dQ, and fix q £ int Q. Let 11+, 11" be two parallel hyperplanes contained 
in R" \ Q and touching dQ from two opposite sides. We also assume that there exists a ball 
B such that Q C B C AB C Uy. Then there exist e^ > (given by Proposition \6.8\) such that, 
if diam((5) < Ec then 

\mr < cinHiQ X y) ""^{d^^^(g>n^)-dist(g,n-)} ^^,,^^^^^,^^^^ 

where £n+ denotes the maximal length among all the segments obtained by intersecting Q with 
a line orthogonal to Il~^ , and 7g (Q ^ ^) is defined as in (23P- 



Proof. Fix q £ Q. Observe that u = on dQ and consider the c-cone h'^ generated by q and 
Q of height —h'^{q) = —u{q) as in (j6.15p . From Lemma l6.7l fd) we have 

\d'h'\{{q})<\d'u\iQ), 

and from Loeper's local to global principle. Corollary 13.41 above, 

dh'{q) = -Dgc{q,d'h'{q)). 

Therefore, 

\dh'\{{q}) < ||detD,VllcO{{,-}xy)|9'^/i1(m). 

The lower bound on |5/i'^|({g}) comes from (|6.16p . This finishes the proof. D 

Combining this with Theorem 13.81 we get the following important estimates: 
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Theorem 6.11 (Alexandrov type upper bound). Assume (B0)-(B3), and let u : Uy i — > R 

be a c-convex function from Theorem \4-^ There exist e^ > small, depending only on the 



dimension and the cost function, and constant C{n), depending only on the dimension, such 
that the following holds: 

Letting Q denote the section {u < 0} C Uy, assume |5'^n| < 1/A in Q and u = on 
dQ. We also assume that there exists a hall B such that Q <Z B C AB C Uy, and that 
diam(Q) < e'^. For ^ < t < 1, let qt ^ Q be a point such that qt G tdQ, where tdQ denotes 
the dilation with the factor t with respect to the center of E, the ellipsoid given by John's 
Lemma (see (|3.9p ). Then, 

\u{qtT < C{n)^{l - t)V2"-^if-(Q)2, (6.24) 

where 7^ = %{Q x V) is defined as in (|4.4p . which satisfies 7^ < J^^^ f'f'om Corollary ^.^ 
Moreover, 

"^"^«^"% C(n)^. (6.25) 



^"(g)2 - ^ ^ A 

Remark 6.12. Aficionados of the Monge-Ampere theory may be less surprised by these 
estimates once it is recognized that the locahzation in coordinates ensures the cost is approx- 
imately aflttne, at least in one of its two variables. However, no matter how well we approxi- 
mate, the non-affine nature of the cost function remains relevant and persistent. Controlling 
the departure from affine is vital to our analysis and requires the new ideas developed above. 

Proof of Theorem 16. Ill For {) < sq < t < 1, lei qt ^ Q he a, point such that qt G tdQ. By 
Theorem 13.81 applied with sq = l/(2n) we can find H"'" 7^ 11" parallel hyperplanes contained 
in T^V \ Q, supporting Q from two opposite sides, and such that 

minjdist (g, n+ ) , dist (g, H" ) } ^ (j(^\n _ i)i/2"-^ 

Then (16.24^ follows from Lemma [6.101 and the assumption \d^u\ < 1/A. 
To prove ()6.25p . observe that, since 2;^Q C g-E'j 



mfQ\( 1^) ur 



< C(n)i(l-— )i/2" 



^"(Q)2 - ' ' A ' 2n 

On the other hand, taking n"*" and 11" orthogonal to one of the longest axes of E and choosing 
q £ 2^ ill Lemma [6.101 yields 

\uiq)r < Cin)^n^^iQf, Vg G ^E. 

A 

Combining these two estimates we obtain (j6.25p to complete the proof. D 
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7 The contact set is either a single point or crosses the domain 

The previous estimates (Theorems 16.21 and I6.1ip may have some independent interest, but 
they also provide key ingredients which we use to deduce injectivity and Holder continuity 
of optimal maps. In this and the subsequent section, we prove the strict c-convexity of the 
c-convex optimal transport potentials u : U i — > R, meaning d'^u{x) should be disjoint from 
d'^u{x) whenever x,x ^ U are distinct. This shows the injectivity of optimal maps, and is 
accomplished in Theorem 18.11 In the present section we show that, if the contact set does 
not consist of a single point, then it extends to the boundary of U. Our method relies on the 
(B3) assumption on the cost c. 

Recall that a point x of a convex set S C R" is exposed if there is a hyperplane supporting 
S exclusively at x. Although the contact set S := d'^ u^ (y) may not be convex, it appears 
convex from y by Corollary 13. 4| meaning its image q{S) C Uy in the coordinates (j4.ip is 
convex. The following theorem shows this convex set is either a singleton, or contains a seg- 
ment which stretches across the domain. We prove it by showing the solution geometry near 
certain exposed points of q{S) inside Uy would be inconsistent with the bounds (Theorem 16.21 
and 16. lip established in the previous section. 

Theorem 7.1 (The contact set is either a single point or crosses the domain). Assume 
(B0)-(B3), and let u be a c-convex solution of (j3.8p with U C U open. Fix x G U and 
y G d'^u{x), and define the contact set S := {x £ U \ u{x) = u{x) — c{x,y)-'rc{x,y)}. Assume 
that S 7^ {x}, i.e. it is not a singleton. Then S intersects dU. 

Proof. To derive a contradiction, we assume that S ^ {x} and S CC U (i.e. S D dU = 0). 

As in Definition 14.11 we transform (x, u) i — > (q, u) with respect to y, i.e. we consider the 
transformation q^Uy i — > x{q) G U, defined on Uy := —Dyc{U,y) C T^V by the relation 

-Dyc{x{q),y) = q, 

and the modified cost function c{q,y) := c{x{q),y) — c{x{q),y) on Uy x V, for which the 
c-convex potential function q £ Uy i — > u{q) := u{x{q)) + c{x{q),y) is level-set convex. We 
observe that c{q, y) = for all q, and moreover the set S = d"^ u^ (y) appears convex from 
y, meaning Sy := —Dyc{S,y) is convex, by Corollarv 13.41 

Our proof is reminiscent of Caffarelli's for the cost c{q, y) = —{q, y) [HI Lemma 3]. Observe 
q := —Dyc{x,y) lies in the interior of the set U^ := —Dyc{U^,y) where l^^^ttl G [A/7^,7^/A], 
according to Corollarv 14.51 Choose the point q^ € Sy CC Uy furthest from q; it is an exposed 
point of Sy. We will see below that the presence of such exposed point gives a contradiction, 
proving the theorem. 

Before we proceed further, note that because of our assumption S CC U, the sets in the 
following argument, which will be sufficiently close to S, are also contained in U; the same 
holds for the corresponding sets in different coordinates. This is to make sure that we can 
perform the analysis using the assumptions on c. 

For a suitable choice of Cartesian coordinates on V we may, without loss of generality, take 
<f ~ Q. parallel to the positive y^ axis. Denote by Cj the associated orthogonal basis for TyV , 
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■DyC{Ko,y) 



-DyCiK^y) 



Figure 7: If the contact set Sy has an exposed point go, we can cut two portions of Sy with two 
hyperplanes orthogonal to q — q^ . The diameter of —Dyc{Ko,y) needs to be sufficiently small to 
apply the Alexandrov estimates from Theorem 16. II [ while ~Dyc{K}^,y) has to intersect U^ in some 
nontrivial set to apply Theorem 16.21 (this is not automatic, as q^ may not be an interior point of 

spt|9^u|). 



and set iP := (g°,ei) and b := (g, ei), so the halfspace {q G T^V ~ R" | gi := (g,ei) > 6°} 
intersects Sy only at g". Use the fact that q^ is an exposed point of Sy to cut a corner Kq 
off the contact set S by choosing s > small enough that h = {1 — s)}fi + sh satisfies: 

(i) —Dyc{Ko, y) := Sy D {q £ Uy \ qi > b} is a. compact convex set in the interior of Uy] 

(ii) diam.{—Dyc{KQ,y)) < s'^/2, where e'^ is from Theorem 16.111 

(iii) —Dyc{KQ,y) C B C 5B C Uy for some ball B as in the assumptions of Theorem 16.111 

Defining q^ := (1 — s)q^ + sq, x^ := x{q'^) the corresponding c-segment with respect to y, 
and q = q^, note that Sy n {qi = b} contains q, and Kq contains x := x^ and x^. Since the 
corner Kq may not intersect the support of \d'^u\ (especially, when q^ is not an interior point 
of sptl^^^ul), we shall need to cut a larger corner Kq as well, defined by —Dyc{KQ,y) := 
Sy n {q £ Uy \ qi > b} , which intersects U^ at x. By tilting the supporting function slightly, 
we shall now define sections K/r C K^ of u whose interiors include the extreme point x^ 
and whose boundaries pass through x and x respectively, but which converge to Kq and Kq 
respectively as e — )■ 0. 

Indeed, set y^ := ?/ + eei and observe 

iTT'lix) := -c(x, ye) + c{x, y) + c(x^ y^) - c(x^ y) 
= e{-Dyc{x, y) + Dyc{x'',y),ei) + o(e) 
= eii-Dycix, y), ei) - (1 - 5)6° - sb) + o{e). (7.1) 
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Figure 8: We cut the graph of u with the two functions m| and ml to obtain two sets K^ ~ Kq 
and K^ « Kq inside which we can apply our Alexandrov estimates to get a contradiction (Theorem 
16. in to K^, and Theorem l6.2l to K^). The idea is that the value of u — ml at xq is comparable to its 
minimum inside K;,, but this is forbidden by our Alexandrov estimates since xq is too close to the 
boundary of Kq. However, to make the argument work we need also to take advantage of the section 
K^, in order to "capture" some positive mass of the c-Monge- Ampere measure. 



Taking s E {s, 1} in this formula and e > shows the sections defined by 

Ke := {x I u{x) < u{x) - c{x, ye) + c(x, Ve)}, 
kI := {x I u{x) < u{x) - c{x, ys) + c(x, ye)}, 

both include a neighbourhood of xq but converge to Kq and Kq respectively as e — > 0. 

We remark that there exist a priori no coordinates in which all sets Ke are simultaneously 
convex. However for each fixed e > 0, we can change coordinates so that both Ke and K^ 
become convex: use ye to make the transformations 

q := -Dyc{xe{q),ye), 



Ce{q,y) ■■= c{Xeiq),y) - c{Xeiq),ye), 



so that the functions 



Ue{q) := u{xe{q)) +c{xe{q),ye) -u{x) - c{x,ye), 
ul{q) := u{xe{q)) + c{xe{q),ye) - u{x) - c{x,ye). 

are level-set convex on Uy^ := Dyc{U,ye). Observe that, in these coordinates, Ke and Ke 
become convex: 

Ke := -Dyc{Ke,ye) = {g G Uy, I Mq) < 0}, 
Kl := -Dyc{Kl,ye) = {(/ G Uy^ \ ul{q) < 0}, 

and either Ke C Ke or Ke C Ke since Ue{q) — ul{q) = const. For e > small, the inclusion 
must be the first of the two since the limits satisfy Kq C Kq and q G Kq \ Kq . 
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In the new coordinates, our original point x G U\ the exposed point x^ , and the c-convex 
combination x with respect to y, correspond to 

qe-=-Dyc{x,ye), q^ := -Dyc{x°,ye), qe := -Dyc{x,ye). 

Thanks to (ii) and (iii), for e sufficiently small we have diam(i^£) < e'^ and K^ C B C AB C 
Uy^ for some ball B (note that this ball can be different from the one in (iii)), so that all the 
estimates of Theorem 16.111 apply. 

Let us observe that, since lim£_5.o Qe — Qe = Q^ ~ Q and q^ E dQ, (j6.24p combines with 
K, C K} and \d^^Ue\iK,) < A-f-^'^iK,) from ^ and Corollary |131 to yield 

l«-<«l" .0 ae.^O. (72) 



A7rJf"(-K'i)2 



On the other hand, x £ S implies Us{q'^) = 'm^{x^), and x £ S implies ul{q'^) = ml{x^) 
similarly. Thus (j7.ip yields 



"Ml?) e(i>» - 6) + o(e) 

Our contradiction with (j7.2p - (|7.3p will be established by bounding the ratio |n^((jr^)|"/^"(ir^)^ 
away from zero. 
Recall that 

5° = (-Z)yc(x°,y),ei) = max{qi \ q £ -Dyc{Ko,y)}) > b 

and u{x) — u{x) > —c{x, y) + c(x, y) with equality at x^. From the convergence of K^ to Kq 
and the asymptotic behaviour (jT.ip of 7nl{x) we get 

uliq'^) -n(x°) - c(x°,y£) + u(x) + c{x,ye) 



inf^i ni sup^g^i[-'u(x(g)) - c{x{q),ye) + ^(x) + c(x,ye)] 

^ -c(x°, i/g) + c(x, ye) + c(x°, y) - c(x, y) 
~ sup^g/^i [-c(x, ye) + c{x, ye) + c{x, y) - c{x, y)] 

^ £{{-Dyc{x'^,y),ei)-b) + o{e) 



(7.4) 



e(max{gi | q G -Dyc{Kj,y)} - b) + o(e) 



1 
> - 

- 2 

for e sufficiently small (because, by our construction, maxjgi | q G —Dyc{Ke,y)} is exactly 
{—Dyc{x^, y), ei)). This shows u^{qe) is comparable to the minimum value of u^. To conclude 
the proof we would like to apply Theorem 16.21 but we need to show that \d'^u\ > A on a stable 
fraction of K^ as e — )• 0. We shall prove this as in |5l Lemma 3]. 

Since K^ converges to Kq for sufficiently small e, observe that K^ (thus, K^) is bounded 
uniformly in e. Therefore the affine transformation (LJ)~^ that sends K^ to Bi C Ke* C Bn 
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as in Lemma 13.71 is an expansion, i.e. \{Ll) ^q — (L^) ^q'\ > Co\q — q'\, with a constant 
Co > independent of e. Since x is an interior point of f/'^, B^g-^ifj {x) C U^ for sufficiently 

small 5 > (here (3~ is from ([O]) ). hence -B25/Co(9e) "^ ^j^,,' where Uy^ := -Dyc{U^,ye). 
Let g* := (L£)~^(g£). Then, by the expansion property of {L\)~^, we have 

C/,V := (Li)-i(t/,^J D B2sm- 

Since iT^' is convex, it contains the convex hull C of -Bi U {q^}- Consider C n B2s{q*)- Since 
dist(0, q*) < n, there exists a ball B* of radius (5/n (centered somewhere in C D B2siqs)) such 
that 



1,* 

e • 



s*ccni?25(gnci?25(g1)ni^, 

Therefore, the ellipsoid 

Es:=Ll{B*) 

is contained in C/r, . Notice that Es is nothing but a dilation and translation of the ellipsoid 
E = LI{Bi) associated to K^ by John's Lemma in (j3.9p . By dilating further if necessary 
(but, with a factor independent of e), one may assume that Eg C B C 4:B C Uy^ for some 
ball B (as before, this ball can be different from that for Kq or -ftTg). Thus, we can apply 
Theorem 16.21 (with Q = K^) and obtain 



inf^iu 



r.lin 



K^'^^^ >^2n 



where the inequality > is independent of e. As e — t- this contradicts (I7.2p - (l7.4p to complete 
the proof. D 

Remark 7.2. As can be easily seen from the proof, when U = U one can actually show 
that if 5 is not a singleton, then Sy has no exposed points in the interior of Uy. Indeed, if by 
contradiction there exists q^ an exposed point of Sy belonging to the interior of Uy, we can 
choose a point q £ Sy sufficiently close to q^ in the interior of Uy = U^ such that the segment 
q^ — q IS orthogonal to a hyperplane supporting Sy at q^. Then it can be immediately checked 
that the above proof (which could even be simplified in this particular case, since one may 
choose Kq = Kq and avoid the last part of the proof) shows that such a point q^ cannot 
exist. 

8 Continuity and injectivity of optimal maps 

The first theorem below combines results of Sections [5] and [7] to deduce strict c-convexity of 
the c-potential for an optimal map, if its target is strongly c-convex. This strict c-convexity 
— which is equivalent to injectivity of the map — will then be combined with an adapta- 
tion of Caffarelli's argument [5] Corollary 1] to obtain interior continuity of the map — or 
equivalently C^-regularity of its c-potential function — for (B3) costs. 
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Theorem 8.1 (Injectivity of optimal maps to a strongly c-convex target). Let c satisfy (BO)- 
(B3) and (B2)s. If u is a c-convex solution of (j3.8p on U^ C U open, then u is strictly 
c-convex on U , meaning d'^u{x) and d'^u{x) are disjoint whenever x, x E U are distinct. 

Proof. Suppose by contradiction that y E d'^u{x) n d'^u{x) for two distinct points x,x E U^, 
and set S = d^ u'^ (y). According to Theorem 17.11 the set S intersects the boundary of U 
at a point x E dU n d'^ u^ {y). Since p.Sp asserts A < \d'^u\ on U and \d'^u\ < A on [/, 
Theorem 1 5. ir a) yields y ^ V (since a;,5; E U ), and hence x E f/ by Theorem 15.1( b). This 
contradicts x E dU and proves the theorem. D 

By adapting Caffarelli's argument [SI Corollary 1], we now show continuity of the optimal 
map. Although in the next section we will actually prove a stronger result (i.e., optimal maps 
to strongly c-convex targets are locally Holder continuous), we prefer to prove this result for 
two reasons: first, the proof is much simpler than the one of Holder continuity. Second, 
although not strictly necessary, knowing in advance that solutions of (j3.8p are C^ will avoid 
some technical issues in the proof of the C^'" regularity. 

Theorem 8.2 (Continuity of optimal maps to strongly c-convex targets). Let c satisfy (BO)- 
(B3) and (B2)s. Ifu is a c-convex solution of (|3.8p on U^ C U open, then u is continuously 
difjerentiable inside U . 

Proof. Recalling that c-convexity implies semiconvexity (see Section [3] and also (14. Sp ). all we 
need to show is that the c-subdifferential d^u{x) of u at every point x ^ U\ \s a, singleton. 
Notice that d^u{x) CC V by Theorem 15.1( a). 

Assume by contradiction that is not. As d^u{x) is compact, one can find a point yo in 
the set d^u{x) such that —Dxc{x,yo) E du{x) is an exposed point of the compact convex 
set du{x). Similarly to Definition 14.11 we transform (x,n) i — > {q,u) with respect to yo, i.e. 
we consider the transformation q E t/yg i — > x{q) E U , defined on Uy^ = —Dyc{U, yo) + 
Dyc{x,yo) C Ty^V by the relation 

-Dyc{x{q),yo) + Dyc{x, yo) = q, 

and the modified cost function c{q,y) := c{x{q),y) — c(x(g),yo) on Uyf-, x V, for which the 
c-convex potential function q E Uy^ i — > u{q) := u{x{q)) — u{x) + c{x{q),yo) — c{x,yo) is 
level-set convex. We observe that c{q, yo) = for all q, the point x is sent to 0, tt > n(0) = 0, 
and u is strictly c-convex thanks to Theorem 18.11 Moreover, since —Dxc{x,yo) E du{x) was 
an exposed point of du{x), = — Dgc(0, yo) is an exposed point of du{0). Hence, we can find 
a vector v E du{0) \{0} such that the hyperplane orthogonal to f is a supporting hyperplane 
for du{0) at 0, see Figure [9j By the semiconvexity of u (see (j4.5p ). this implies that 

u{-tv) = o(t) for t > 0, u{q) > {v, q) - Mc\q\'^ for ah q € Uy^. (8.1) 

Let us now consider the (convex) section K^ := {u < e}. Recalling that the unique 
minimum point with u(0) = 0, we have that K^ shrinks to as e — >• 0, and u > on 
Uyg \ {0}. Thus by (18. ip it is easily seen that for e sufficiently small the following hold: 

KeC{q\ {q, v) < 2e}, -a{e)v E Ke, 
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Figure 9: w e du{0) and the hyperplaiie orthogonal to v is supporting du{0) at 0. 




s/o(l) 



Ib 



Figure 10: Since tire hyperplane orthogonal to v is supporting du{0) at 0, we have u{—tv) = o{t) for 
t > 0. Moreover, by the semiconvexity, u grows at least linearly in the direction of v. 



where a(e) > is a positive constant depending on e and such that Q;(e)/e — )• +oo as e — )■ 0. 
Since is the minimum point of tt, this immediately implies that one between our Alexandrov 
estimates (|6.4p or (j6.24p must be violated by u inside Kg, for e sufficiently small, which is the 
desired contradiction, see Figure [101 D 



9 Engulfing property and Holder continuity of optimal maps 

Our goal is to prove the C^^ regularity of u inside U . The proof will rely on the strict 
c-convexity of n inside U (Theorem El| and the Alexandrov type estimates from Theorems 
16.21 and 16.111 Indeed, these results will enable us to show the engulfing property under the 
(B3) condition (Theorem 19. 3p . thus extending the result of Gutierrez and Huang |29l Theo- 
rem 2.2] given for the classical Monge- Ampere equation. Then, the engulfing property allows 
us to apply the method of Forzani and Maldonado |^ to obtain local C^'" estimates (The- 
orem [93|). Finally, by a covering argument, we show that the Holder exponent a depends 
only on A and n, and not on particular cost function (Corollary 19. 6p . 



Remark 9.1. We point out that in this section the assumptions (B2) and (B2)s are used 
only to ensure the strict c-convexity (Theorems 18. ip . Since the following arguments are 
performed locally, i.e. after restricting to small neighborhoods, having already obtained the 
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strict c-convexity of u, the c-convexity (B2) and (B2)s of the ambient domains are not 
further required. This is a useful remark for the covering argument in Corohary 19.61 



Given points x G U^ and y G d'^u{x), and r > 0, we denote by S{x,y,T) the section 

S{x, y, t) := {x £ U \ u{x) < u{x) — c{x, y) + c{x, y) + r}. (9-1) 

Notice that by the strict c-convexity of u (Theorem 18. ip . S{x,y,T) — >• {x} as r — )• 0. 

In the following, we assume that all points x, x that we choose inside U are close to each 
other, and "relatively far" from the boundary of U\ i.e. 

dist(2;, x) < min ( dist(x, dU^),disi{x, dU^)) . 

This assumption ensures all the relevant sets, i.e. sections around x or x, stay strictly inside 
U\ 

As we already did many times in the previous section, given a point y G d'^u{x) with 
X G U\ we consider the transformation (x, n) i — > {q,u) with respect to y (see Definition 
14.11). and define the sections 



Qr = {g G Ty*V I u{q) < r}, r > 0. 

Note that Q^ corresponds to S{x, y, r) under the coordinate change. By the (B3) condition, 
each Qr is a convex set (Theorem 14. 3|) . We also keep using the notation pQr to denote the 
dilation of Q^ by a factor p > with respect to the center of the ellipsoid given by John's 
Lemma (see (13. 9p ). 

In our analysis, we will only need to consider the case r ^ 1. This is useful since, thanks 
to the strict c-convexity of u (Theorem 18. ip . we can consider in the sequel only sections 
contained inside U"^ and sufficiently small. In particular, for every section Q = S{x, y, r) or 
Q = S{x,y,KT), after the transformation in Definition 14.11 

- we can apply Theorem 16.21 with Eg = E (so (|6.4p holds with 5 = 1); 

- we can apply Lemma [6.31 with J^ = Q; 

- Theorem 16.111 holds . 

The following result generalizes |291 Theorem 2.1(ii)] for c{x,y) = —x ■ y to (B3) costs: 

Lemma 9.2 (Comparison of sections with different heights). Assume (B0)-(B3) and let u 

be a strictly c-convex solution to (13. 8p on U . Take x G U , y G d'^u{x) and r sufficiently 
small so that 

Sr=S{x,y,T)cU^ 

and set Qr := —Dyc{Sr,y)- Then, there exist < pQ < 1, depending only on the dimension 
n and ^^Ic /^ (^*^ particular independent of x,y and t), such that 

Qt/2 ^ PoQt- 
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Proof. This is a simple consequence of Theorems 16.21 and 16. Ill indeed, considering u — t, for 
p > l/2n ([OIjl gives 



while by ([O]) (with (5 = 1) 
Hence, we get 



if"(QO' < C(n)^T". 
A 



\uiq)-T\<C{n)^{l-p)' 



<C{n) 



-^ _ 2 

A 



(1 - P? 



where the last inequality follows from jfj^ /A^ < [7(^7c /A] as in Corollary 14.51 Therefore, 
for 1 — po sufficiently small (depending only on the dimension n and ^^Ic /^) ^^ S^* 

\u{q) -t\ <t/2 yq€Qr\ poQr, 

so Qr/2 ^ PoQt as desired. D 

Thanks to Lemma [9.21 and Lemma [6. 3[ we can prove the engulfing property under (B3), 
extending the classical Monge- Ampere case of |29( Theorem 2.2]: 

Theorem 9.3 (Engulfing). Assume (B0)-(B3) and let u be a strictly c-convex solution to 
([3l8]) on U^. Let x,x £U^ be close: i.e. 



dist(x, x) <C min ( dist(x, dU^),dist{x, dU^)) . 

Then, there exists a constant K > 1, depending only on the dimension n and 7^7~/A, such 
that, for all y G d'^u{x), y G d'^u{x), and r > small (so that the relevant sets are in U^), 

x £ S{x,y,T) =^ X £ S{x,y,KT). (9.2) 

Proof. First, fix {x,y) G d'^u. We consider the transformation {x,u) i — > {q,u) with respect 
to y (see Definition 14. ip . Let q= —Dxc{x,y) G T^V denote the point corresponding to x in 
these new coordinates. To show (19.21) we will find K > such that 



q£Q^ =^ u{q)<uiq) + c{q,y)-c{q,y)+KT, Vy G 9S(g). (9.3) 

Fix T > small, q G Qt and y G d'^u{q). Assume by translation that the John ellipsoid of 
Q2t is centered at the origin. By Lemma [9. 2 1 we have Q^ C P0Q2T for some po < 1 (depending 
only on the dimension n and ^flc /^)- Hence, by Lemma 16.31 applied with J^ = Q2T we 
obtain 

\\Dqciq,y)\\*Q^^<c(n,^^)T VgGQ., (9.4) 
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where || • ||q is the dual norm associated to Q2t (see (j6.5p ). 

Recah from Theorem 14.31 that q minimizes u, so u{q) < u{q). To obtain (I9.3P from this, 
we need only estimate the difference between the two costs. Since both q,q & Qr C Q2t, by 
the definition of || • ||q and (19.4p . (recalling sq + {1 — s)q G Qr due to convexity of Qr) we 
find 

\c{q,y) -c{q,y)\ 

/ Dgc{sq+ {1- s)q,y)ds- {q- q) 
Jo 

/ Dqc{sq+{1- s)q,y) ■ qds + / Dgc{sq + {1 - s)q,y) ■ q 
Jo Jo 

<2 I \\Dgc{sq+{l-s)q,y)\\*Q^^ds 
Jo 



< 



ds 



<C7(n, ^"^" 



A 
thus establishing ()9.3p . 



D 



Having established the engulfing property, the C^'" estimates of the potential functions 
follows by applying a modified version of Forzani and Maldonado's method [24) . Here is a 
key consequence of the engulfing property. 



Lemma 9.4 (Gain in c-monotonicity due to engulfing). Assume the engulfing property (j9.2 
holds. Let c,u,x,x and K he as in Theorem \ 9. 31 and let y £ d^u{x) and y G d'^u{x). Then, 



l + K 
K 



u{x) - u{x) - c{x, y) + c(x, y)] < c(x, y) - c(x, y) - c(x, y) + c(x, y). 



Proof. Given x, x, notice that u[x) — u{x) + c{x, y) — c(x, y) > 0. Fix e > small, to ensure 

r := u{x) — u{x) + c(x, y) — c(x, y) + e > 0. 
Then x € S{x,y,T), which by the engulfing property implies x € S{x,y,KT), that is 

u{x) < u{x) + c(x, y) — c{x, y) + K[u{x) — u{x) + c{x, y) — c(x, y) + e]. 
Letting e — )• and rearranging terms we get 

{K + l)u{x) <{K + l)n(x) + c(x, y) - c(x, y) + K[c{x, y) - c(x, y)], 
or equivalently 



u{x) — u{x) < 



l + K 



[c{x,y) -c{x,y)] + 



K 

l + K 



[c{x,y) -c{x,y)]. 
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This gives 

K 

u{x) - u{x) - c{x, y) + c{x, y) < ^ [c(x, y) - c(x, y) - c{x, y) + c(x, y)], 

as desired D 

In the above lemma, it is crucial to have a factor (1 + K)/K > 1. Indeed, the above 
result implies the desired Holder continuity of n, with a Holder exponent independent of the 
particular choice of c (see Corollary 19. 6p : 

Theorem 9.5 (Holder continuity of optimal maps to strongly c-convex targets). Let c satisfy 
(B0)-(B3) and (B2)s. If u is a c-convex solution of (|3.8p on U C U open, then u G 
^loc (^ )' with K as in Theorem 1 9. 51 which depends only on the dimension n and j^j" /X. 

Proof. As we already pointed out in the previous section, although not strictly needed, we 
will use the additional information that u £ C^{U'^) (Theorem 18. 2p to avoid some technical 
issues in the following proof. However, it is interesting to point out the argument below works 
with minor modifications even if u is not C^, replacing the gradient by subdifferentials (recall 
that u is semiconvex), and using that semiconvex function are Lipschitz and so differentiable 
a.e. We leave the details to the interested reader. 

The proof uses the idea of Forzani and Maldonado |24| . The c-convexity of u is strict on 
U\ according to Theorem 18.11 Given a point Xs € U\ we denote by ys the unique element 
in d'^u{xs)', the uniqueness of ys follows from the C^ regularity of u, since yg is uniquely 
identified by the relation Du{xs) = —Dxc{xs,ys)- 

Let xq G U . We will show that for xi E U sufficiently close to xq, 



\u{xo) - u{xi) - Du{xi) ■ {xq - Xi)\ < \xo - Xl 



l+l/K 



from which the local C^'^' regularity of u follows by standard arguments. 
Fix a direction v with \v\ small, set Xg = xq + sv, and consider the function 

(p{s) := u{xs) - u{xo)+c{xs,yo) - c(xo,yo)> 0, 

for s £ [0, 1]. The idea is to use Lemma [9.41 to derive a differential inequality, which controls 
the growth of (p. First, observe that 

(j)'{s) = Du{xs) ■ V + Dxc{xs,yQ) ■ v. 

Since Du{xs) = -Dxc{xs,ys), we get 

(l)\s)s = [D^c{xs,yo) - Dxc{xs, ys)] ■ {sv) 

> c{xo,ys) -c{xs,ys) - c{xo,yo) + c{xs,yo) - \\Dl^c\\Lac(_u^Y)s'^\v\'^ . 

So, by Lemma [9.41 we get 

^ 4>{s) < (p'is)s+ \\Dl^c\\L-^(uxV)s'^\vf , 
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that is 

d ( 0(s) "^ \\DIxC\\l^{UxV)\'v\ 



dt Vsi+V^y - s^/K 

Hence (^(sj/s^+i/^^ < (/.(I) + \\DI^c\\lo.(uxV)\v? j] ^'^'^ dr < 0(1) + Ci (since 1 - ^ > 0). 
So, 

(l){s)=u{xs) - u{xo) - c{xo, yo) + c{xs,yo) 

< s^-^^/'^[u{xi) - uixo) - c(xi) + c{xo + V, yo) + Ci] 
< 2CiS "*" ' (choosing \v\ small enough and using the continuity of u and c). 

By the arbitrariness of Xo,v and s we easily deduce that, for all xq, xi G U sufficiently close, 

u{xi) -u{xo) -c{xo,yo) + c{xi,yo) < 2Ci\xo -xil^+^Z^^. 

Since c is smooth, Du{xi) = —Dxc{xi, yi), and u{xq) — u{xi) — c(xi, yo) + c(xo, yo) ^ 0, the 
last inequality implies 

\u{xo) - u{xi) - Du{xi) ■ {xq - Xl)\ 

< \u{xo) - u{xi) -c{xi,yo) + c(xo,yo)| + \\Dl^c\\Lco(^uy,v)\xo - 2;iP 

<C72(|xo-xi|i+i/^ + |xo-xi|2) 

<2C2|xo-xi|^+^/^ 

for all Xo,xi G U sufficiently close. This proves the desired estimate, and concludes the 
proof of the C^^J regularity of u inside U . D 

In fact, the Holder exponent in the previous theorem does not depend on the particular 
cost function: 

Corollary 9.6 (Universal Holder exponent). With the same notation and assumptions as in 
Theorem \9.5\. u G ^^^{U^), where the Holder exponent a > depends only on n and A > 0. 

Proof. Recalling 7^ = j^{U-'^xV) := ||(det D^j^c)^-^||^oo((7Axy) from (I4.4J) . we see 7+ — > 1/7" 
as the set U x V shrinks to a point. Since u is C^ by Theorem 19. 5( d'^u gives a single- valued 
continuous map. Compactness of U combined with (BO)-(Bl) allows the set d'^ur\{U x V) 
to be covered with finitely many neighborhoods C/fcX V^, such that 7(^(f/feX T4)7~(C/fcX V^) < 2 
for all k. Hence, thanks to Remark 19.11 we can apply Theorem l9.5l on each such neighborhood 
([/^ n f/fc), which in turn yields a Holder exponent < a < 1 depending only on n and 
A >0. D 
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